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TNO  I  Kennis  voor  zaken 

Responsie  van  verschillende  munitieartikelen 
bij  opwarming  en  brand 


Brand  geeft  op  zichzelf  al  gevaarlijke  siluaties.  Als  die  in  combinatie  met 
munitie  optreedt.  worden  de  gevaren  nog  groter.  De  resultaten  van  dit 
rapport  geven  de  mogelijkheid  de  kritische  condi ties  van  munitie  bij  brand 
te  bepalen  ter  ondersteuning  van  het  brandbestrijdingproees. 


Probleems  telling 

Veiligheid  van  munitie  voor  bet  eigen 
personeel  en  mate  rice  l  onder  a  lie  omstan- 
digheden  is  een  van  de  voorwaarden  voor  het 
liinctioneren  van  de  krijgsmachi.  Naar  aan- 
leiding  van  de  vuurwerkramp  in  Enschede 
wordl  er  nu  binnen  de  Kk  gekeken  naar  de 
veiligheid  van  munitieopslagen.  Het  doel  van 
dit  project  is  bet  vaststellen  van  het  gedrag 
van  munitie anikelen  in  brand,  Het  lesultaat 
kan  als  basis  diene n  vtxir  aanhevelingen  voor 
een  her/iening  van  de  riehllijn  MP40  21  en 
h  rand  weemple  idi  ng . 

Beschrijving  van  de  werkzaamheden 

Het  project  is  go  start  met  het  oniwikkden 


van  een  rekenmodel  om  de  cook-off  tijden 
van  munitieartikelen  te  bepalen,  Met  dit 
mkenmodel  kan  een  sehatting  gemaakt 
worden  van  de  lijd  en  temperatuur  tot 
expkrsie  en  bet  zogenaarnde  'point-of-no- 
retum’ .  Van  eendrietal  munitieartikelen  zijn 
deze  lijd  en  temperatuur  tot  explosie  bepaald 

Resultaten  en  concltisies 

[ >e  tijden  en  temperaturen  tot  explosie  van  de 
volgende  munitie  anikelen  z.ijn  berekend: 

35  mm,  1 20  mm  en  de  1 55  mm,  Vervolgens 
is  een  verge  lij  king  gemaakl  van  betekende 
waaiden  met  literatuurwaarden  betretTende 
thermische  response  van  an  ike  Jen, 

LX'  kritische  temperatuur  waarbij  een 


munilieamke!  zal  exploderen  hangl  sterk  at 
van  de  opwarmsnelheid,  Voor  de  mcesie 
conventionele  verschoten  munitie  met  de 
expkrsieven  HMX.  RDX  en  TNT  als  hoold- 
bestanddeel  en  een  diameter  in  de  range  van 
35-155  mm  ligt  deze  grens  op  een  lijn  tussen 
1 50  bij  een  opwarming  van  0,001  TVs 
(3.6  °CAiur)  naar  200  bij  een  opwarming 
van  04  TVs  (360  °C/uur)  en  opwarm- 
snelheden  daar  boven,  Ite  lijd  die  daarmee 
gepaard  gaat  is  aJhanketijk  van  de 
temperatuur  op  moment  van  atlezen  en  de 
opwarmsneiheid,  Kruiten  rcageren  sneller. 
Voor  artikden  met  kniit  ligt  deze  grens  iager, 
van  105  °C  naar  150  °C  bij  sndheden  van 
0.00 1  K/stot  1  K/s,  Voor  munitie  in  een 
direct e  brand  is  uiterste  voorzichligheid 
gebodeiL  Tijden  van  minder  dan  1  minuet 
zijn  niet  vrcemd.  Om  een  beter  overzicht  van 
de  kruiten  te  krijgen  zal  een  gmtere  variatie 
van  lsnel  rcagetende  kruiten1  bekeken 
moeten  worden,  De  incest  gangbare 
cxplosieven  in  conventionele  munitie  (RDX. 
HMX,  LIST  en  Comp  B  tRDX/TNT),  Comp 
A3  (RDX)  etc.)  kunnert  op  de  zojuisl 
genoemde  manier  benaderd  worden.  Dit 
geldt  wk  v*x>r  vee!  modeme  PBXen  met  als 
hoofdbestanddeel  RDX.  HMX  of  TNT  en 
bindersystemen  als  LTITB  Ook  bijvtxrrrieeld 
de  modeme  155  munitie  van  Rheinmeiall 
(met  KS32)  vah  binnen  deze  normen  net  als 
munitie  met  het  veelgebruikte  blast*- 
expkrsief  PBX-N 109. 
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bckond  /ijn.  I  )ii  kan  gemote  n  worden  met 
bchulp  van  dckmimstfic  tempo  rat  nur¬ 
se  nsoren  in  do  bunkers.  Verier  meet  bekend 
zijn,  hoeveel  on  welke  soorten  munitie 
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I  Inleiding 


Ongelukken  in  het  recente  verleden,  zoals  de  vuurwerkramp  in  Enschede,  hebben  de 
aandacht  gericht  op  de  opslag  van  explosieve  stoffen  en  muni  Heart  ikelen  en  de  gevaren 
die  dit  met  zich  meebrengt.  Door  deze  ongelukken  zijn  de  verantwoordelijke  personen 
opnieuw  gewezen  op  deze  gevaren.  Specifiek  bij  een  brand  in  een  munitiemagazijn 
dient  de  brandweer  bij  het  verkennen  en  het  b lessen  van  de  brand  veilig  ingezet  te 
kunnen  worden.  Daarvoor  is  het  noodzakelijk  dat  zij  bekend  zijn  met  het  gedrag  van  de 
opgeslagen  muni  tie  in  een  brand. 

De  richtlijnen  voor  munitiebrandbestrijdtng  zijn  vastgelegd  in  MP  40-21,  'Voorschrift 
voor  opslag  gevaarlijke  stoffen’  en  de  Brandweeropleiding  Brandmeester  [5,  6]. 

De  laatste  keer  dal  deze  richtlijnen  zijn  aangepast  is  na  een  aantal  studies  in  de  periode 
van  1993-1995  [7,  8],  In  veiligheids studies  wordt  aange nomen  dat  bij  een  brand  in  een 
munitiemagazijn  alle  aanwezige  muni  tie  betrokken  is  en  maximaal  zal  reageren/ 
detoneren,  Er  wordt  dus  van  een  'Maximum  Credible  Event'  {MCE}  ultgegaan. 

Dit  betekent  dat  maximale  en  dus  conservative  veiligheidsafstanden  in  acht  genomen 
worden.  In  het  laatste  decennium  is  de  kennis  op  het  gebied  van  munitie  en  explosieve 
stoffen,  en  speciaal  het  gedrag  bij  een  ongewenste  reactie,  sterk  toegenomen.  Veel  van 
deze  ontwikkel ingen  komen  voon  uit  het  onderzoek  naar  Insensitive  Munitions 
(Minder  Kwetsbare  Munitie,  MKM).  Het  doel  van  MKM  is  dat  de  munitie  minder 
heftig  reageen  als  zij  ongewenst  tot  reactie  komt.  Deze  ontwikkelingen  betekenen  dat 
veel  meer  informatie  gegeven  kan  wwden  over  de  heftigheid  van  de  reactie  van  munitie 
in  een  brand  en  de  tijdsperiode  die  verstrijkt  voordat  deze  reactie  optreedt.  Tevens  is  het 
tegenwoordig  niet  meer  noodzakelijk  om  in  alle  gevallen  van  een  MCE  uit  te  gaan. 

Deze  informatie  is  van  groot  belang  voor  de  verkenning  en  bestrijding  van  een  muni  tie- 
brand.  Om  deze  redenen  heeft  het  Hoofd  Brandwreerzorg  van  het  NATCO  te  Gouda  de 
behoefte  geuit  om  in  kaart  te  brengen  wat  nu  het  reete  brandgedrag  en  -effect  van 
(een  bulk)  munitie  kan  zijn  wanneer  deze  betrokken  zijn  bij  een  brand  (cook-off). 

Binnen  dit  project  is  voor  een  aantal  verschillende  mu  ni  heart  ikelen  in  opslag, 
afhanketijk  van  het  soon  op  warming  en  daaraati  gekoppeld  de  op  warms  nelbeid  van 
deze  artikelen,  de  temperatuur  tot  reactie  bij  een  brand  numeriek  bepaald.  Om  de 
respons  van  de  munitieartikelen  te  bepalen  zijn  eerst  fysische  eigenschappen  van  de 
stoffen  in  de  munitie  in  de  literatuur  opgezocht.  Voor  zover  er  literatuurgegevens  over 
de  respons ie  van  de  munitie  zijn,  zijn  deze  verge leken  met  de  bepaalde  waarden. 

Op  deze  manier  is  er  een  grens  bepaald  waarbij  de  munitie  we  I  of  n  iet  reageert 

Hoofdstuk  2  gee  ft  een  samenvatting  van  de  bezoeken  aan  de  Afdeling  Beproeving  van 
Wapensystemen  en  Munitie  (ABWM)  in  Oldebroek  en  het  MunitieMagazijnComplex 
(MMC)  in  Veenhuizen  om  een  idee  te  krijgen  welke  munitie  op  welke  wijze  opgeslagen 
ligt.  Hoofdstuk  3  beschrijft  het  model  en  de  parameters  die  gebruikt  zijn  om  de  tijd  en 
temperatuur  tot  explosie  te  bepalen,  Hoofdstuk  4  geeft  de  resultaten  van  de 
berekeningen  weer  van  de  geseleeteerde  munitie,  Hoofdstuk  5  toont  de  resultaten  van 
een  gevoeligheidstudie  met  betrekking  tot  de  invoerpara meters.  Hoofdstuk  6  kijkl  naar 
de  verschiilen  van  twee  RDX  houdende  explosieve  stoffen  waarvan  den  met  64  wt % 
RDX  en  een  met  85  wt%  RDX,  in  vergelijking  tot  expert mentele  gegevens  van  cook-off 
readies,  Hoofdstuk  7  kijkt  naar  de  gevaren  bij  bulkopslag  van  met  name  kruiten. 
Hoofdstuk  8  besluit  met  conclusies  en  aanbevelingen. 


TNO-rapport  TNO-DV  2007  A315 


6/41 


2  Bezoeken  munitiemagazijnen 


()m  een  goed  bee  Id  te  krijgen  van  hoe  muni  tie  opgeslagen  ligt,  zijn  de  Afdeling 
Beproeving  van  Wapensystemen  en  Munitie  (ABWM)  in  Oldebroek  en  het 
MunitieMagazijnComplex  (MMC)  in  Veenhuizen  bezocht.  Een  beknopte  samenvatting 
van  deze  bezoeken  staan  in  de  volgende  paragrafen,  een  uitgebreider  verslag  is  te 
vinden  in  referentie  [2].  De  manier  waarop  met  brand  omgegaan  wordt,  kan  per  land 
verschillen  aangezien  elk  land  zelf  verantwoordelijk  is  voor  bet  eigen  brandweerbeleid. 

2.1  Afdeling  Beproeving  van  Wapensystemen  en  Munitie 

Bij  de  afdeling  Beproeving  van  Wapensystemen  en  Munitie  worden  wapensystemen  en 
munitie  beproefd  op  hun  functioneren.  Er  worden  onder  andere  fuel  fine  proeven, 
classificatietesten  periodiek-  en  calamiteitenonderzoek  uitgevoerd. 

De  voorraden  munitie  aanwezig  op  hei  complex  in  Oldebroek  liggen  opgeslagen  in 
bunkers,  zie  figuur  2.  L  Per  bunker  staat  op  de  buitenkant  welke  klasse  er  opgeslagen 
ligt  (klasse  1:  massaexpiosief,  etc.).  De  huidige  regels  over  banrieres  en  afstanden  die  op 
het  terrein  gelden  zijn  in  de  MP40-21  richttijnen  vastgelegd.  De  bunkers  zijn  over  het 
algemeen  genomen  niet  vol:  de  munitie  is  in  kleine  hoeveelheden  aanwezig.  Wei  zijn  er 
veel  verschillende  soorten  munitie  aanwezig,  zie  figuur  2,2,  De  munitie  die  tijdens  de 
testdag  gebruikt  wordt.  wordt  in  kleinere  hoeveelheden  opgeslagen  als  dagopslag  in  een 
aparte  ruimte  in  het  be wer kings gebouw,  De  hoeveelheden  munitie  die  opgeslagen  ligt, 
is  bekend  en  wordt  bijgehouden 


Figuur  2, 1  Opslagbunker  bij  de  afdeling  Beproeving  van  Wapensystemen  en  Munitie  in  Oldebroek. 
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Ftguur  2,2  De  bunker  ligi  diet  geheel  vol,  versehillende  test  mate  rialen  worden  hier  opgeslagen. 

In  geval  van  een  brand  melding  wordt  het  terrein  overgedragen  aan  de  gemeentelijke 
brandweer,  De  gemeentelijke  brandweer  is  veramwoordelijk  voor  het  blussen  van  de 
brand.  De  brandweer  wordt  ingelicht  door  iemand  van  het  ABWM  over  de  plaatsen 
waar  welke  soon  munitie  ligt.  Op  die  moment  zijn  de  afstanden  die  gehouden  rnoeten 
worden  in  gevat  van  brand  bepaald  op  basis  van  de  vergunningen  over  hoevee)  er 
maximaal  mag  liggen  en  niet  op  basis  van  de  hoeveelheden  die  er  uerkelijk  aanwezig 
zijn.  In  geval  van  brand  ontmimen  de  mensen  van  het  ABWM,  De  brandweer  gaat  wel 
het  terrein  op,  maar  houdt  zich  aan  de  geldende  richtlijnen;  de  procedure  in  geval  van 
brand  in  een  magazijn  is  1 000  m  afstand  houden  en  wachten. 

Bij  de  Afdeling  Beproeving  van  Wapensystemen  en  Munitie  brengt  vooral  het 
bewerken  en  het  testen  van  de  munitie  risico's  met  zich  mee.  Wanneer  de  munitie 
opgeslagen  ligt  in  een  hunker  is  er  maar  een  kleine  kans  dat  er  iets  gebeurt.  Bij  het 
ABWM  is  er  gesproken  orn  in  elke  ruirnte  camera's  aan  te  brengen,  De  videobeelden 
kunnen  dan  laten  zten  wat  er  gaande  is  in  de  ruimtes.  Helaas  bleek  dat  dit  systeem  veel 
te  duur  was  (ongeveer  EUR  6,000  000,-)  om  te  kunnen  uitvoeren. 

2.2  MunitieMagazijnComplex 

In  Nederland  zijn  er  op  verschillende  plaatsen  munitiemagazijneomplexen.  De  munitie 
ligt  opgeslagen  in  bunkers,  De  belading  van  de  bunkers  is  bekend  via  het  MUNIS 
systeem,  Wanneer  de  veranderingen  van  een  bunker  groter  dan  10%  zijn*  wordt  het 
systeem  bijgewerkt.  In  de  bunkers  in  Veenhuizen  ligt  meer  opgeslagen  dan  in  de 
bunkers  in  Oldebroek,  Figuur  2.3  laat  een  voile  opslagbunker  zien. 
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In  elke  bunker  hangt  een  overzieht  wat  er  niaxiniaal  aan  gewieht  van  een  bepaalde 
klasse  mag  liggen.  Daarnaast  bevinden  zich  in  elke  bunker  de  procedures  voor 
ontruiming  en  de  procedure  wat  te  doen  bij  brand.  De  verpakking  waarin  de  muni  tie 
opgeslagen  ligt.  is  niet  voor  elk  artikel  gelijk:  bet  meest  gebruikt  is  hout  geVmpregneerd 
met  brandvertragers. 


Figuur  2.3  In  do  bunkers  in  Veenhuizen  ligt  meer  rnunitie  opgeslagen  dan  in  do  bunkers  in  Oklebirak. 

Voor  transport  wordt  de  rnunitie  in  containers  ge laden.  De  bunkers  zijn  zo  groot  dat 
beladen  containers  in  de  bunker  passen,  zie  figuur  2.4.  De  containers  worden  maar  met 
1  laag  rnunitie  geladen.  Om  te  voorkomen  dat  de  rnunitie  gaat  schuiven  wordt  een 
container  helemaal  ingcbouwd  met  hout.  Op  deze  tranier  kan  de  container  ook  op  zijn 
kop  geplaatst  worden  zonder  dat  er  schade  ontstaat  aan  de  artikelen,  Op  de  plants  van 
bestemming  wwdt  de  rnunitie  uit  de  container  gehaald  en  opgeslagen  volgens  de 
richtlijnen. 
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Figuur  2.4  Containers  met  muniiie  die  klaar  staan  voor  transport  worden  tx>k  opgeslagen  in  bunkers. 

De  brand  weer  op  bet  complex  heeft  een  aanvalsplan  geschreven  over  hoe  te  handelen 
bij  brand.  Hierin  staat  vermeld  dat  er  gekeken  moet  worden  waar  de  betreffende  brand 
is  en  welke  klasse  muni  tie  er  opgeslagen  ligt.  Dit  aanvalsplan  is  zowel  bij  de  bedrijfs- 
brandw^eer  als  bij  de  gemeentelijke  brand  weer  bekend.  Een  aantal  mensen  van  de 
bedrijfsbrandweer  zit  ook  in  de  gemeentelijke  brandweer,  hierdoor  is  bet  contact  tussen 
beide  brandweerdiensten  goed.  Het  aanvalsplan  hevat  veiligheidscirkels  die  in  acht 
genomen  worden  bij  brand.  Ook  hier  geldt  de  MP40-21 

Op  het  terrein  in  Veenhuizen  zijn  twee  werkplaatsen.  Wanneer  er  muni  tie  aanwezig  is, 
mo  gen  de  werkplaatsen  alleen  betreden  worden  met  speciale  (anti-statisehe)  schoenen. 
Bij  elke  werkplek  is  een  brandblusser  aanwezig.  De  muren  tussen  de  werkplekken  zijn 
zo  ontworpen  dat  deze  moeten  blijven  staan  in  geval  van  een  explosie.  De  achterwand 
wordt  wel  naar  buiten  geblazen.  Bij  iedere  werkplek  hangt  een  camera  om  te  registreren 
wat  er  gebeurt. 

2.3  Conclusies  munitieopslagplaatscn 

In  Nederland  wordt  op  verschillende  plaatsen  muni  tie  opgeslagen.  De  bunkers  op  de 
opslaglocaties  liggen  aanzienlijk  voller  dan  de  bunkers  bij  afdelingen  waar  onderzoek 
gedaan  wordt,  De  richtlijnen  van  de  MP40-21  gelden  op  alle  locaties.  Tijdens  vredes- 
operaties  worden  de  regels  in  acht  gehouden  zoals  beschreven  in  de  materieellogistieke 
instructie  "Munitieveiligheid  tijdens  vredesoperaties\ 

Op  dit  moment  zijn  er  in  de  bunkers  geen  camera’s  en  temperatuursensoren,  die  op 
afstand  uitgelezen  worden,  Het  wordt  aanbevolen  om  wel  camera's  en  temperatuur- 
sensoren  te  plaatsen  zodat  op  afstand  mogelijke  op  warming  geidentificeerd  kan  worden 
en  de  aard  en  locatie  van  deze  opwarming  bepaald  kan  worden,  Na  aflezing  van  de 
geregistreerde  temperatuur  en  de  daaruit  afgeleide  opwarmsnelheid  kan  eventueel 
besloten  worden  of  er  nog  aktie  genomen  wrordt  om  de  brand  te  blussen  of  dat  er  op  een 
afstand  van  1QOQ  m  gewacht  wordt  tot  dat  een  eventuele  explosie  heeft  plaatsgevonden. 


TNO-rapport  |  TNO-DV  2007  A315 


10/41 


3  Gebruikt  model  en  parameters 


Qrn  de  responsietemperatuur  van  verschillende  munitieartikelen  te  bepaien  wordt  een 
l-dimensionaal  model  gebruikt,  Dtt  model  heeft  verschillende  eigenschappen  van  de 
stoffen  in  de  muni  tie  nodig.  In  paragraaf  3. 1  wordt  het  model  kort  besproken. 

In  paragraaf  3.2  warden  de  parameters  van  de  stoffen  besproken. 

3,1  Gebruikt  mode) 

Wanneer  een  munitieartikel  in  of  nabij  een  brand  ligt,  zal  deze  opwarmen.  Door  deze 
opw^arming  kan  het  energetische  materiaal  in  het  artikel  gaan  reageren.  Wanneer  er  door 
deze  reactie  meer  warmte  geproduceerd  wordt  dan  dat  er  aan  de  omgeving  afgestaan 
wordt.  zal  het  artikel  opwarmen.  ontsteken  en  eventueel  exploderen,  Deze  opwarming 
kan  op  verschillende  manieren  gesimuleerd  worden.  In  Scholtes  [3]  en  [4]  worden 
verschillende  modellen  beschreven  waaronder  het  I-dimensionale  simulatiemodel, 
geVmplementeerd  in  het  spreadsheetprogramma  Microsoft  Excel  2002.  De  modellen  zijn 
ontwikkeld  voor  zowel  de  Koninklijke  Landmacht  als  voor  de  KomnkUjke  Marine. 

Het  l-dimensionale  model  bevat  4  verschillende  soorten  berekeningen,  namelijk4 
verschillende  manieren  van  op  warming,  De  eerste  is  opwarming  door  middel  van  vrije 
convectie;  ofwel  een  natuurlijke  opwarming  aan  de  lucht.  De  tweede  manier  van 
opwarming  is  op  warming  door  middel  van  direct  contact,  bijvoorbeeld  een 
verwarmingsspiraal  om  het  artikel  been.  De  derde  manier  is  de  bepaling  van  het  'point 
of  no  return*.  Hierbij  wordt  de  mini  male  temperatuur  gezocht  waarbij  het  artikel  nadat 
de  opwarming  van  buitenaf  gestopt  is  (bijvoorbeeld  de  brand  is  geblust)  toch  nog  zijn 
ontsteektemperatuur  bereikt.  De  laatste  vorm  is  opwarming  door  middel  van  straling 
van  vuur.  Figuur  3.1  geeft  de  lay-out  van  de  spreadsheet  weer. 

ten  ander  gebruikte  methode  is  een  bestaand  l-dimensionaal  pakket  dat  door  TNO 
ontwikkeld  is  (GSCO)  waarmee  direct  opwarming  van  munitie  doorgerekend  kan 
worden,  Helaas  is  het  met  deze  code  niet  mogelijk  een  raket motor  door  te  rekenen  of 
opw' arming  aan  de  lucht  die  langzaam  warmer  wordt 
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Figuur  3- 1  V oorbecld  v  an  de  Exce  1  -spreadsheet * 


De  spreadsheet  heeft  de  stofeigenschappen  van  het  energetische  materiaal  en  reken- 
parameters  nodig.  In  label  3. 1  staan  de  rekenpara meters.  De  stofeigenschappen  worden 
besproken  in  paragraaf  3,2, 


Tabet  3.1  Rekenparameters  voor  hoi  model 


Parameter 

Symbooi 

Eenheid 

Waarde 

Opwarmsnelhefd 

P 

K/s 

varieert:  1, 0*1,  0,01 
en  0.001 

W  arm  teo  v  e  rd  rac  htsc  o§ff  i  c  i£n  t 

h 

W/(m2K) 

8 

Mogelijkheid  tot  opwarming 

optie  1 

- 

1  of  0 

Opwarming  via  vrije  convectie  of 
direct  contact 

optie  2 

- 

Oof  1 

Mogelijkheid  voor  point  of  no  return 

optie  3 

- 

Oof  1 

Mogelijkheid  voor  de  stralingswarmte 

optie  4 

- 

0  of  1 

Straal  van  mogelijke  Hbore,r 

Rb 

m 

afh.  v.  artikel 

Straal  van  explosief' 

Rc 

m 

afh.  v.  artikel 

Straal  van  liner 1 

Ri 

m 

afh,  v.  artikel 

Tijdstap 

t 

s 

Varieert 

Opstarttijd  van  het  vuur 

tstart  lot  vuur 

s 

60 

Startle  mperatuur 

T**t 

K 

293*15 

Omgevingstemperatuur 

Tang 

K 

293,15 

Temperatuur  van  het  vuur 

Tt 

K 

1300 

Emissiecoefficient  van  het  vuur 

£r 

- 

0,9 

Emissiecoefficient  van  het  oppervlak 

£5 

- 

0,8 

i 


In  de  bo  re  ken  ingen  is  aangenomen  dat  do  ci  Under  uir  altcen  RDX  bestaai  en  gecn  ‘bore'  of  liner  be  vat 
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Tijdens  de  berekeningen  wordt  gezocht  naar  het  tijdstip  waarbij  bet  explosief  zijn 
ontsteektemperatuur  bereikt.  Wanneer  ergens  in  bet  explosief  de  ontsteektemperatuur 
van  dll  explosief  bereikt  wordt,  wordt  de  wandtemperatuur/oventemperatuur  afgelezen. 
Bij  de  directe  opwarming  wordt  de  wandtemperatuur  opgelegd  en  bij  op  warming  via 
convectie  wordt  de  omgevingstemperatuur/oventemperatuur  opgelegd.  Zie  figuur  3.2 
voor  een  illustratie  van  de  wand-  en  oventemperatuur.  De  opwrarming  via  convectie  laat 
de  op  warming  van  reele  situaties  zien. 


I'iguur  3,2  Op  warming  via  direct  contact  =  opgdegde  wandtemperatuur  (links)  en  via  convectie  - 
opgelegde  oragevingstemperatuur  (rechts). 


3.2  Stofafhankelijke  eigenschappen  ais  input  voor  mode] 

3.2.1  Exp  lost  even 

In  het  model  zijn  eigenschappen  van  de  stof  nodig.  Voor  explosieven  zijn  veel  van  deze 
gegevens  bekend.  In  tabel  3.2  staan  de  gebruikte  waarden  voor  de  stoffen  die 
voorkomen  in  granaten,  waarvan  de  tijd  tot  explosie  en  de  temperatuur  van  de  wand 
berekend  is. 


Tabel  3.2  Stofafhankelijke  eigenschappen  van  de  explosieven  gebruikt  in  de  berekeningen. 


Parameter 

Symbool 

Eenheid 

TNT 

Ref-  [11] 

RDX  (RU81) 
Ref.  [12] 

HMX  (HU27) 

Ref.  [12] 

Warmtegeleid¬ 

baarheid 

A 

W/(m-K) 

0,2092 

0,25 

0,317 

Dichtheid 

P 

kg/m3 

1570 

1580 

1640 

Warmtecapaciteit 

Cp 

J/(kg-K) 

1500 

1360 

1277 

Reactiewarmte 

Q 

J/kg 

1255200 

2088890-0, 85 1 

2343400  0,87' 

Pre-exponentlSIe 

coefficient 

2 

1/s 

2,510-10" 

3,16228-1 0ts 

5t0-1019 

Active  rings- 
energie 

Ea 

J/mol 

143900 

199000 

220290 

Ontsteek¬ 

temperatuur 

TofllSJSftk 

aC 

297 

203 

275 

;  De  composite  bestaat  voor  deze  Fractie  uii  het  explosief. 


3.2.2  Kruiten 

Om  de  responsietemperatuur  en  tijd  te  bepalen  van  de  voortdrijvende  lading  en  de  klein 
kaliber  munitie  is  het  nodig  om  de  stofeigenschappen  van  het  kruit  te  weten. 
Eigenschappen  zoals  de  warmtegeleidbaarheid.  k.  en  de  diehtheid,  p,  zijn  hierbij  lastig 
te  bepalen.  De  warmtegeleidbaarheid  wordt  namelijk  beinvloed  doorde  lucht  die  tussen 
het  kruit  zit.  De  vorm  van  het  kruit  heeft  hier  ook  invloed  op,  Wanneer  het  kruit  dichter 
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op  elkaar  gestort  is,  zal  de  warmtegeleidbaarheid  meer  en  meer  lijken  op  de  waarde  van 
de  warmtegeleidbaarheid  van  bet  kruit  zelf  Om  een  goed  beeld  te  krijgen  moet  in  het 
model  gerekend  worden  met  een  waarde  die  ligt  tussen  die  van  massief  kruit  en  die  van 
lucht,  Een  aannemelijk  range  voor  de  geleidbaarheid  is  van  0X>6  tot  0T3  W/mK. 

Voor  de  dichtheid  geldt  een  gelijke  redenatie.  De  lucht  tussen  het  kruit  zorgt  ervoor  dal 
de  schijnbare  dichtheid  van  het  artikel  lager  is  dan  de  dichtheid  van  het  kruit  zelf. 
Aannemelijke  waarde  n  voor  de  dichtheid  zijn  200  tot  900  kg/m\  De  warmtecapaciteit 
varieert  iets  minder  en  Ligt  in  de  range  van  1300  tot  1500  i/kg  K, 

De  parameters  zoals  de  pre-exponentiele  factor  en  de  activeringsenergie  zijn  iastiger 
vast  te  leggen  in  een  range.  De  pre-exponentiele  factor  en  de  activeringsenergie  zijn  aan 
elkaar  gekoppeld  en  mogen  niet  zomaar  los  van  elkaar  gevarieerd  warden, 

Een  bij  ko  mend  aspect  is  dat  kruit  (nitrocellulose)  op  zichzelf  geen  stable le  stof  is, 

De  -O-NO:  binding  breekt  waardoor  er  radicalen  vrijkomen,  Deze  radicalen  reageren 
vervolgens  weer  met  een  ander  nitro-oxide  groep  waardoor  uiteindelijk  het  kruit 
thermisch  wegloopt.  Om  deze  reactie  te  vertra gen/ voo rkome n  wordt  er  een  stabilisator 
toegevoegd  aan  het  kruit,  De  hoeveelheid  stabilisator  wordt  in  de  loop  van  de  tijd 
verbmikt  De  stabilisator  wordt  verbruikt  door  het  afvangen  van  de  radicalen  en 
voorkomt  op  deze  manier  de  kettingreactie.  Door  het  toevoegen  van  een  stabilisator 
wordt  de  activeringsenergie  van  het  kruit  verhoogd,  Voor  de  bereken ingen  is  het 
belangrijk  om  te  weten  hoe  de  activeringsenergie  (gekoppeld  aan  de  pre-exponentiele 
factor)  zich  verhouden,  Wat  nu  de  waarde n  van  het  kruit  zijn.  is  moeilijk  te  bepalen, 
Deze  worden  beinvloed  door  de  stabilisator,  De  waarden  van  de  stabilisator  zijn  dus  van 
groot  belang  omdat  ze  de  reactiekinetiek  veranderen.  In  de  Uteratuur  zijn  weinig 
gegevens  bekend  over  de  reactiekinetiek  van  kruiten  al  dan  niet  met  stabilisator. 

Daarom  is  het  lastig  om  'juiste*  waarden  te  bepalen  voor  de  berekening  van  de 
responsie  van  hel  kruit,  Hierbij  komt  nog  dat  de  hoeveelheid  stabilisator  die  het  kruit 
bevat  op  het  moment  van  cook-off  niet  bekend  is  en  in  de  praktijk  natuurlijk  varieert. 

De  eigensc happen  van  het  kruit  kunnen  varieren  in  de  orde  van  groottes  zoals  vermeld 
in  label  3,3.  In  overleg  met  de  toenmalige  projectbegeleider  ir  M  Pot  is  besloten  om  in 
dit  rapport  een  handbreedte  te  bepalen  door  verse  hi  llende  eigensc  happen  in  te  schatten. 
Door  een  'worst  case’  en  een  "best  case'  te  maken,  kan  een  brandbreedtc  van  de 
responsie  gemaakt  worden.  Voor  de  activeringsenergie,  de  reactie  war  mte  en  de  pre- 
exponentiele  factor  zijn  waarden  geschal  met  behulp  van  de  gemiddelde  warmte 
produetie  van  kruiten  gemeten  in  Heat  Flow  Calorimetry  (HFC)  experimenten  bij 
85  °C.  Deze  warmteproductie  is  gecombineerd  met  een  activeringsenergie  wpat  leidt  tot 
4  situaties.  Men  zou  verwachten  dat  bij  een  hoge  activeringsenergie  het  kruit  het  meest 
stabiel  is.  Echter  uit  de  berekeningen  bleek  dat  een  hoge  activeringsenergie 
gecombineerd  met  een  lage  warmteproductie  bij  85  °C  leidde  tot  een  hogere 
energieproductie  dan  bij  een  lage  re  activeringsenergie,  Daarom  zijn  beide  cases 
doorgerekend.  Bij  de  hoge  warrnteproductie  bij  85  DC  geldt  hetzelfde,  ook  hier  zijn 
beide  cases  berekend.  De  overige  eigenschappen  van  kruit  best  case  L  en  2  zijn  a  Is  een 
heel  gunstig  kruit  genomem  uit  het  oogpunt  voor  de  cook-off  Dit  zijn  de  hoogste 
waarden  voor  de  warmtegeleidbaarheid.  de  dichtheid  en  de  warmtecapaciteit.  Voor  de 
eigenschappen  van  kruit  worst  case  1  en  2  zijn  hiervoor  juist  de  laagste  waarden 
genomea 
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Tabel  3.3  Slofeigenschappen  van  kruiten. 


Symbool 

Eenheid 

Best  case  1 

Best  case  2 

Worst  case  1 

Worst  case  2 

Warmtegeteidbaarheid 

h 

W/rrvK 

0.3 

0.3 

0,06 

0,06 

Dichtheid 

P 

kg/m3 

900 

900 

200 

200 

Warmtecapaciteit 

Cp 

J/kg-K 

1500 

1500 

1300 

1300 

Pre-exponenti£te 
coefficient  Reactiewarmte 

ZQ* 

W/kg 

5,518  1014 

2,443-1 01fl 

4,415-1  O'* 

1 .954  IQ19 

Activeringsenergle 

Ea 

J/mot 

110000 

135000 

110000 

135000 

Ontsteektemperatuur 

4  antsteek 

•c 

160 

160 

160 

160 

De  waarden  van  Z  en  Q  ztjn  samcngenomcn. 


Deze  geschatte  waarden  zijn  verge le ken  met  de  waarden  uit  de  literatuur.  Tabel  3.4  laat 
de  eigenschappen  van  kruit  KB6954  en  kruit  C  zien  die  overgenomen  zijn  uit  [  1  ]. 

Kniit  M  bevat  de  eigenschappen  die  overgenomen  zijn  uit  [10]. 


Tabel  3,4  Stofeigemehappen  van  kruiten. 


Parameter 

Symbool 

Eenheid 

KB6954 

[1] 

Kruit  C 

m 

Kruit  M 
[10] 

Warmtegeleidbaarheid 

A 

W/m  K 

0,08 

0,08 

0.08 

Dichtheid 

P 

kg/m3 

840 

840 

840 

Warmtecapaciteit 

Cp 

J/kg-K 

1350 

1350 

1350 

Pre-exponentiSle 

co§ffici6nt-Reactiewarmte 

Z  Qr 

W/kg 

5,782- 10t6 

1,608  IQ'7 

3,54- 1018 

Activeringsenergie 

Ea 

J/mot 

122700 

120000 

188200 

Ontsteektemperatuur 

Toni  sleek 

°C 

160 

160 

160 

J  De  waarden  van  Z  en  Q  zijn  same n genome n. 


Over  het  cook-off  gedrag  van  kruiten  is  weinig  kennis  beschikbaar.  Het  schatten  of 
bepalen  van  de  parameters  die  nodig  zijn  om  het  gedrag  te  voorspellen  is  lastig. 

Echter  over  de  levensduur  van  kruiten  is  meer  bekend,  Om  een  voorspelling  over  de 
levensduur  te  doen  word!  er  gekeken  naar  de  warmteproductie  bij  85  °C  zoais  hierboven 
vermeld.  In  combi natie  hiermee  wordt  er  ook  gekeken  naar  het  verbruik  van  de 
stabilisator  Om  een  schatting  van  het  cook-off  gedrag  te  kunnen  doen.  kan  ook  bepaald 
worden  hoe  lang  het  duurt  tot  de  stabilisator  verbruikt  is.  Aangenomen  wordt  dal  het 
kruit  gaat  reageren  wanneer  de  stabilisator  verbruikt  is.  Om  de  levensduur  te  bepalen 
wordt  de  tijd  gesehat  aan  de  hand  van  vergelijking  (1 ). 


[-E. 

! 

--T 

l  * 

Js 

Tn  J  j 

waarin; 

t*  =  veilige  levensduur  bij  omgevingstemperaluur  [$]; 
tm  =  duur  van  de  meting  [s]; 

T*  =  o mge vingS'/ops  1  a g te mperatuur  [K]; 

Tm  =  meettemperatiiur  [K]. 

Hierbij  wordt  verder  aangenomen  dat  de  warmteproductie  en  de  stabilisatorconsumptie 
verlopen  volgens  de  wet  van  Arrhenius,  het  kruit  de  temperatuur  van  de  omgeving  heeft 
en  dat  het  react iepad  bij  vcrsehillende  temperaturen  gelijk  is. 
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De  voortdrijvende  lading  wordt  geinitieerd  door  eon  slaghoedje  met  een  vlatnpijp. 

De  vlampijp  bevat  voornamelijk  nitrocellulose,  de  grootste  component  in  het  kruit. 
De  compositie  in  het  slaghoedje  is  over  het  algemeen  een  initiaalspringstof, 

De  springstoffen  zijn  gevoelig  voor  stoot  en  soms  ook  wrijving.  Echter  op  thermisch 
gebied  zijn  deze  stable  1.  Hierdoor  zal  door  thermische  stimuli  het  mitiatiegedeeke 
waarschijnlijk  pas  later  ontsteken  dan  het  kruit.  Daarom  zijn  deze  stoffen  buiten  de 
beschouwing  van  dit  project  gehouden. 
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4  Berekeningen  aan  geselecteerde  munitie 

Er  worden  veel  verschillende  artikelen  gebruikt  bij  de  landmacht.  Er  is  een  selectie  van 
de  meest  gebruikte  munitie  gemaakt  om  de  hoeveelheid  berekeningen  beperkt  te 
houden.  De  volgende  artikelen  zijn  geselecteerd:  de  Howitzer  munitie  (155  mm),  de 
1 20  mm  mortier  en  de  35  mm  patronen.  Zowel  de  Howitzer  als  de  1 20  mm  bevatten  een 
granaat  met  een  explosieve  stof  alswel  een  voortdrijvende  lading  bestaande  uii  kruit. 

Bij  de  Howitzer  is  dit  kruit  verpakt  in  kardoezen,  terwijl  bij  de  120  mm  patronen  bet 
kruit  in  een  huls  zit  en  bij  de  120  mm  mortieren  bestaat  bet  kruit  uit  ringen. 

De  berekeningen  aan  de  granaten  van  de  Howitzer  en  120  mm  mortier  worden  in 
paragraaf  4. 1  besproken.  Paragraaf  4.2  bespreekt  de  berekeningen  aan  bet  kruit. 

4.1  Berekeningen  aan  explosieven 

4. LI  Granaat  van  155  mm 

De  1 55  mm  munitie  (Howitzer  munitie)  bestaat  uit  een  granaat  met  TNT  als  explosieve 
stof,  die  weggeschoten  wordt  met  behulp  van  een  kardoes  Het  kardoes  be  vat  kruit, 

De  thermische  responsie  van  het  kruit  wordt  besproken  in  paragraaf  4.2. 


Piguur  4.1  155  mm  granaten  worden  verschoten  mei  een  Howitzer. 
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[-iguur4.2  155  mm  granaai. 

Otn  de  explosietijd  en  temperatuur  van  de  granaai  te  berekenem  is  de  granaai 
gesimuleerd  als  een  cilinder  bestaande  uit  TNT  met  en  zonder  een  sialen  omhulling. 

De  stofeigenschappen  van  TNT  zijn  overgenomen  van  Victor  [11].  In  label  4,  l  staan  de 
gegevens  voor  de  verschillende  stralen,  De  cilinder  wordt  gesimuleerd  als  een  massieve 
TNT  cilinder  daarom  is  de  straal  van  de  'bore'  (gebruikt  bij  berekeningen  aan  een  raket) 
gelijk  aan  0.  De  straal  tussen  de  propellant  en  de  liner  is  zo  verdeeld  dat  de  ruinate- 
stappen  in  de  berekening  gelijk  verdeeld  worden.  Tijdens  de  berekeningen  wordt 
gezocht  naar  het  tijdstip  waarbij  het  expiosief.  in  dit  geval  TNT,  zijn  ontsteek- 
temperatuur  bereikt. 


label  4,  l  Invoergegevens  voor  dc  straal  voor  de  155  mm  granaai  met  TNT, 


Parameter 

Eenheid 

Waarde 

R* 

m 

0 

RP 

m 

0,0568 

R, 

m 

0,0775 

Tabel  4.2  laat  expiosietemperaturen  zien  bij  verschillende  simulaties.  de  opwarming 
vindt  plaats  via  convectie  of  direct  contact  (conv  versus  direct)  met  of  zonder  sialen 
wand  (sw  versus  kaal).  Naast  het  gebruikte  model  zoals  heschreven  in  hoofdstuk  3  is  er 
een  tweede  model  (GSCO)  gebruikt.  Dit  model  (1-dimensionale  TNO  code  geschreven 
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in  Delphi/Pascal)  is  beperkt  (kan  alleen  op  warming  via  direct  contact*  maar  be  vat  wel 
een  stale n  wand),  maar  is  ter  verificatie  van  de  resultaten  van  het  ontwikkelde  model 
opgenonien.  De  berekening  met  de  opwarming  via  convectie  zonder  wand  (kaal)  gaf 
geen  betrouwbare  resultaten  en  is  daarom  weggelaten  in  het  rapport*  De  wand  die  hij 
echte  muni  tie  aanwezig  is  heeft  te  veel  invloed  op  de  berekeningen,  De  Tno  return  is  de 
temperatuur  waarbij  het  artikel  nadat  het  uit  bet  vuur  gehaald  is  toch  nog  tot  reactie 
komt.  De  temperatuur  van  het  point  of  no  return  is  moeilijk  te  bepalen  met  de 
ontwikkelde  spreadsheet.  Ook  bier  heeft  de  wand  van  staal  een  grote  invloed  op  deze 
berekening.  Met  de  Ddimensionale  code  GSCO  is  deze  bepaald  via  directe  opwarming, 
De  temperature n  van  het  point  of  no  return  liggen  echter  dicht  bij  de  temperatuur 
waarbij  het  artikel  reageert.  Bij  langzamere  op  warming  liggen  deze  temperaturen  iets 
lager  niaar  het  is  in  de  praktijk  niet  aan  te  raden  om  nog  in  te  grijpen  wanneer  het  artikel 
deze  temperatuur  bereikt  heeft,  Echter  het  bewustzijn  dat  het  artikel  nadat  het  deze 
temperatuur  heeft  bereikt  alsnog  kan  exploderen  terwijl  de  opwarming  van  buitenaf  al 
wfeg  is,  is  belangrijk,  fn  paragraaf  4.1,2  zal  dit  aan  de  hand  van  figure n  geillustreerd 
worden. 


Tabel  4.2  Temperatuur  van  de  wand  op  moment  van  explosie  bi  j  verschillende  opwarmsnelheden 
top  warming  door  convectie  (conv)  en  direct  contact  (direct)  met  en  /under  stalen  wand 
(sw  versus  kaal)  en  met  behulp  van  GSCO)  voorde  155  mm  granaat  met  TNT. 


3  [K/s] 

"l"sw  pcmv  rci 

7 sw  direct  [°0] 

T kaal  direct  [  C] 

Tosco  sw  direct  PC) 

T GSCO  no  rolurn  PC] 

1 

298 

447 

432 

301 

301 

0r1 

248 

270 

270 

264 

259 

0,01 

211 

216 

220 

219 

217 

0,001 

173 

177 

177 

175 

172 

In  figuur  4.3  is  een  voorbeeld  te  zien  van  het  temperatuurverloop  van  de  155  mm 
granaat  bij  een  op  warming  van  0.1  K/s  door  vrije  convectie*  waarbij  de  temperatuur  Ti 
de  temperatuur  in  de  buitenste  gesimuleerde  cel  is  en  Tm.|o  de  temperatuur  van  de 
binnenste  gesimuleerde  cel  in  de  granaat*  Tijdens  de  berekening  wordt  er  geen  rekening 
gehouden  met  het  smelten  en  eventueel  koken  van  het  TNT, 


TowffrVwill 

- cak:  TwbII 

T| 

—  TM 

- T|-2 

- - Tf-S 

- Tm 

Tm-i 

- Tnv? 

- Tm-3 

Tm-4 
TnvJ 
Tm6 
—  Tm-7 

Tnv9 

- Tnvfl 

Tiri'fO 


Figuur  4.3  Voorbeeld  van  een  resultaat  van  een  berekening  aan  een  1 55  mm  granaal  gevuld  met  TNT 
waarhij  hel  artikel  dtvor  luchl  opgewarmd  word!  (OJTYs). 
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Figuur  4  4  geeft  de  e xp los ie te mperature n  als  functie  van  de  opwarmsnelheid  grafisch 
weer.  Het  is  duidelijk  te  zien  dat  de  methode  via  de  ‘staler*  wand  berekening'  (Tsw  con) 
de  laagste  waarden  voor  de  verschillende  opwarmsnelheden  op  te  vert-  De  waarden 
liggen  allemaal  boven  deze  donkerblauwe  curve  (voorzien  van  een  ruit).  Dit  is  de  worst 
case- 


500 


too 

50 

o 

0,001  0,01  oj  1 

Heating  rate  [K/s] 

i'iguur  4,4  Explosive mperatuur  van  de  155  mm  grarmat  met  TNT  als  tunctie  van  de  opwarmsnelheid  aan 
de  hand  van  de  verschillende  methode  n. 


Voor  extreem  snelle  op  warming  (artikel  direct  in  vuur)  is  het  beter  om  naar 
berekeningen  te  kijken  via  straling  van  het  vuur.  In  tabei  4.3  is  hiervan  het  resultant 
gegeven.  Hierbij  is  straling  de  belangrijkste  term  en  loopl  binnen  60  seconden  de 
temperatuur  op  naar  1300  K.  Hierdoor  zal  het  artikel  naar  schatting  in  71  seconden 
reageren.  Echter  hierbij  en  ook  bij  bovenstaande  berekeningen  moet  nog  in  acht 
genomen  worden  dat  TNT  bij  80  graden  smelt.  Dit  is  met  in  de  berekeningen 
meegenomen  (te  complex),  Dit  kan  eventueel  nog  vertragend  werken  maar  zal  niet 
extreem  zijn. 


Tabei  4.3  Tijd  L‘n  temperatuur  tot  explosie  hij  op  warming  via  directe  warmtestraling  van  vuur.  voor  de 
1 55  mm  granaat  met  TNT, 


PfK/sl 

t&KpSoafl  [$] 

rci 

Twand  [°C] 

N.v.t. 

71 

1027 

434 

4.  L2  Granaat  van  de  120  mm 

De  120  mm  munitie  bestaat  globaal  uit  twee  delen  met  een  aparte  lading;  een  gedeelte 
met  een  springlading  en  een  gedeelte  met  een  voortdrijvende  lading.  Figuur  4.5  en  4.7 
laten  een  120  mm  mortier  zien,  de  ringen  met  kruit  zijn  duidelijk  te  7,ien.  Figuur  4.6  laat 
een  aantal  120  mm  patronen  zien.  Het  bovenste  deel  is  de  granaat  en  het  onderste  deel  is 
de  huls  met  de  voordrijvende  lading,  De  voortdrijvende  lading  wordt  in  paragraaf  4.2 
besproken.  De  springlading  komt  in  deze  paragraaf  aan  bod.  Als  voorbeeld  wordt  de 
lading  van  een  DM12  patroon  genomen,  De  lading  bestaat  uit  een  mengsel  van  1,8  kg 
Compositie  B  (RDX  en  TNT).  Het  patroon  heeft  een  dikwandig  stalen  granaatlichaam 


TNO-rapport  |  TNO-DV  2007  A315 


20/41 


Figuur  4,5  Mortier  1 20  mm. 


-liAir.T.  &UlJ-S£» 


^mt(W33Al 


TJM,  DM1* 


Figuur4.6  Voorbeelden  van  120  mm  patrcmen. 


Figuur  4.7  Voorbeetd  van  een  1 20  mm  mortier. 

Frota  en  Bailey  [13]  hebben  thermisehe  analyses  gedaan  met  mengsels  van  RDX  en 
TNT.  De  resultaten  geven  aan  dat  mengsels  in  de  verhoudingen  40:60  tot 
0:100  RDX:TNT  het  gedrag  vertonen  van  TNT,  Mengsel  in  de  overige  verhoudingen 
(100:0  tot  40:60)  vertonen  het  gedrag  van  RDX.  Composite  B*  de  for  mule  ring  van  de 
granaat  in  een  1 20  mm  patroon  bestaat  uil  60%  RDX  en  40%  TNT.  Aan  de  hand  van 
bovenstaande  reference  kan  aangenomen  worden  dat  composiiie  B  gesimuleerd  kan 
warden  aisof  het  alleen  uit  RDX  bestaat.  Om  een  vergelijking  te  maken  tussen  RDX  en 
TNT  is  de  responsie  voor  TNT  ook  bepaald. 
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label  4.4 

Lnvoergegevens  voor  de  1 20  mm  granaat  mot  RDX  of  TNI, 

Parameter 

Eenheid 

Waarde 

Rb 

m 

0 

RP 

m 

0,044 

R, 

rrt 

0,060 

Tabel  4.5  laat  explosietemperaturen  zien  bij  verschillende  simulaties.  Dc  opwarming 
vindt  plaats  via  convectie  of  direct  contact.  Net  als  bij  de  vorige  berekening  is  bet 
tweede  model  (GSCO)  gebruikt  ter  verifies tie. 


Tabel  4.5  Temperatuur  van  de  wand  op  moment  van  explosie  bij  verschillende  opwarmsne Ibeden 

(opwannmg  door  convectie  (conv)  en  direct  contact  (direct)  met  en  /under  Staten  wand  (swt  en 
met  behulp  van  GSCO)  voor  dc  120  mm  granaai  met  RDX, 


3 

[K/S] 

Tsw  conv 

(“Cl 

T  sw  direci 

[»C] 

Tdjrflci 

rc] 

T QSCO  sw  direel 

[°CI 

Tno  reium  scuaarfcfcrwt 

[°C] 

T-,0  raurn  GSCO 

rc] 

1 

231 

310 

302 

256 

242 

242 

0.1 

204 

216 

217 

217 

218 

214 

0,01 

178 

183 

188 

186 

202 

184 

0,001 

159 

158 

159 

158 

157 

158 

In  figuur  4.8  zijn  de  temperaturen  van  de  explosie  voor  een  120  mm  granaat  voor  de 
verschillende  methoden  bij  elkaar  gezet,  De  waarden  Itggen  allemaal  boven  de  blauwe 
curve  (  voorzien  van  een  ruit).  Dit  is  de  worst  case. 


350 


100 


50 

0 

0,001  DfQ1  0t1  1 

Heating  rate  [K/s] 

Iriguur  4.&  hxplosietemperatuur  van  dc  1 20  mm  granaat  met  KDX  als  t'unctie  van  dc  opwarmsnelhcid  aan 
do  hand  van  do  verschillende  methoden. 

Voor  extreem  snelle  opwrarming  (artikel  direct  in  vuur)  is  het  beter  om  naar 
berekeningen  te  kijken  via  straling  van  het  vuur.  In  tabel  4.6  is  hiervan  het  resultant 
gegeven.  Hierbij  is  straling  de  belangrijkste  term  en  loopt  binnen  60  seconden  de 
temperatuur  op  naar  1300  K,  Hierdoor  zat  het  artikel  naar  schatting  in  56  seconden 
re  age  re  n. 
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Tabet  4.6  Tijd  cn  temperatuur  tot  explosie  via  opwarming  door  straling  van  een  brand,  voofdc  120  mm 
gran  a  at  inoi  RDX, 


3  [K/s] 

t*xp]osiKi  [s] 

rc] 

Twjmd  [°C] 

N.v.t 

57 

976 

291 

De  Tntire[Ljm  is  de  temperatuur  waarbij  het  artikel  nadat  het  uil  hei  vuur  gehaald  is  nog 
net  explodeert.  Met  de  waarden  nit  label  4.5  zijn  figuur  4.9,  4.10  en  4.1 1  gemaakt. 
Figuur  4.9  laat  het  temperatuurprofLel  zien  waarbij  het  art  ike  l  nog  net  niet  reageert  bij 
een  opwarming  van  0,1  K/s,  dit  is  1  °C  onderTnoretun,.  Figuur  4. 10  laat  het  profiel  zien 
waarbij  het  artikel  net  explodeert  bij  een  opwarming  van  0,1  K/s.  Bij  de  opwarm- 
snelheid  van  0,001  K/s  is  na  opwarming  tot  een  temperatuur  van  157  °C  ofwel  430  K 
(zie  gratiek  44  1  en  tabel  4.5),  de  temperatuur  van  het  explosief  vrijwel  homogeen. 
Doordat  het  artikel  uit  de  warmte  gehaald  wordt,  koelt  de  buitenwand  al  af  en  zakt  deze 
temperatuur,  terwijl  binnen  in  het  explosief  door  zel  fop  warming  het  explosief  steeds 
warmer  wordt  tot  dat  een  explosie  optreedt,  zie  figuur  4. 1 1 ,  De  waarden  van  GSCO  en 
de  spreadsheetberekening  liggen  hier  vrij  goed  bij  elkaar,  behalve  bij  de  waarde  van 
0.01  *C/s+  Deze  waarde  is  te  hoog  uitgevallen  in  de  spreadsheet  omdat  er  meer 
rekenstappen  nodig  zijn  in  een  spreadsheet  (>16000,  dus  rekenkundig  probleem  van  de 
spreadsheet).  Bovenstaande  figuren  iliustreren  dat  het  goed  is  om  bewust  te  zijn  van  het 
feit  dat  het  artikel,  nadat  het  deze  temperatuur  heeft  bereikt,  a  Is  nog  kan  exploderen 
terwijl  de  opwarming  van  buitenaf  al  weg  is. 


Tijd  [s] 


TqvwVwbII 

—  -etk:  TwiUi 

Tl 

• — Tl-1 

- Tl-2 

—  Tl-3 

—  Tm 
Tffiht 

- Tm-? 

- Tm-3 

Tm’* 
Tm-5 
Tm-B 
—  Trrv? 
Ttri-8 

—  Tm.ft 
Tm-10 


Figuur  4ry  Het  artikel  koelt  weer  af.  {3  —  0,1  K7s,  opwarming  via  convert  ie  voor  de  1 2u  mm  granaat  met 
RDX. 
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Tijd  [s] 


Tovwrt'waul 

— -calc  Turn 
T1 

- TM 

- Thfl 

- Tl-3 

- Tm 

Tfli-t 
— -  Tm-? 
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tm-fl 
- — -Tm-9 
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Figuur  4  JO  Nad  at  hci  ariikel  nit  hot  vuur  is*  reageert  hot  ulsnog.  p  =  QJ  K/s,  op  warming  via  convectie  voor 
do  1 21 J  mm  granaat  met  RDX. 
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Rguur  4, 1 1  Nadat  hei  artikel  uk  het  vuur  is,  reageert  het  abnog,  p  =  0,001  K/s.  terwijl  de  wand  aan  hot 
afkoelen  is;  opwarming  via  convectie  voor  de  120  mm  granaai  met  Rl>X 


Berekeningen  met  TNT 

Tabel  4.7  laat  explosietemperaturen  zien  bij  verse  hi  llende  si  mu  lattes,  de  op  warming 
vindt  plaats  via  convectie  of  direct  contact.  Net  als  bij  de  vorige  berekeningen  is  het 
tweede  model  (GSCO)  gebruikt  ter  verificatte. 
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l  abel  4,7  Temperatuur  van  de  wand  op  moment  van  explore  bij  verechillende  opwarmsnelhedeD 

(Opw arming  door  eonveetie  (conv)  en  direct  contact  (direct)  met  en  zonder  stalen  wand  (sw>  en 
met  behulp  van  GSCO)  voorde  120  mm  granaat  met  TNT 


p 

m 

T  sw  conv 

[°C] 

"few  evx:t 

[X] 

"Tkaal 

[X] 

T GSCO  sw  direct 

[X] 

T  ,o  return  spreadsheet 

[X] 

T.o  return  GSCO 

[°C] 

1 

288 

394 

386 

300 

296 

300 

0.1 

247 

261 

262 

262 

265 

258 

0.01 

204 

213 

219 

215 

240 

211 

0,001 

177 

175 

175 

175 

173 

173 

In  figuur  4.12  is  de  temperatuur  van  de  explosie  voor  een  1 20  mm  granaat  met  TNT 
voor  de  verschillende  methoden  bij  elkaar  gezet,  De  waarden  liggen  aflemaal  boven  de 
blauwe  curve  (voorzien  van  een  ruit )*  Dit  is  de  worst  case. 


450 


100 

50 

0 

0,001  0,01  0,1  1 

Heating  rate  [aC/s] 

Figuur  4, 12  Hxptosieternperaturen  als  timelie  van  de  opwarmsneiheid  voor  TNT  voor  een  1 >  mm  granaat. 

De  T^rtfum  temperature n  zijn  in  tabel  4*7  gegeven  en  verge leken  met  de  res u bate n  van 
GSCO.  De  resultaten  van  deze  laatste  liggen  dicht  bij  de  explosietemperaturen. 

Ook  hier  geldt  dezelfde  conclusie  als  voor  de  RDX  berekeningen.  De  point  of  no  return 
temperaturen  liggen  erg  dicht  bij  de  reaetietemperaturen,  Men  zal  reken ing  moeten 
houden  met  bet  alsnog  reageren  van  de  artikelen  als  deze  wel  tot  de  point  of  no  return 
temporal  uur  opgewarmd  zijn  ter  w  i j  t  daarna  de  op  warming  a!  gestopt  is, 

Voor  extreem  s nolle  opwarming  (artikel  direct  in  vuur)  is  het  nog  beter  om  naar 
berekeningen  te  kijken  via  straling  van  het  vuur  In  tabel  4.8  is  hiervan  het  resultaat 
gegeven.  Hierbij  is  duidelijk  dat  al  na  SI  seconden  een  explosie  kan  optreden,  Echter 
hierbij  en  ook  bij  bovenstaande  berekeningen  moet  nog  in  acht  genome n  worden  dat 
TNT  bij  80  graden  smelt.  Dit  is  met  in  de  berekeningen  meege nomen  (te  complex). 

Dit  kan  eventueel  nog  vertragend  werken.  maar  zal  niet  extreem  zijn. 


Tabel  4.8  Tijd  en  temperatuur  tot  explosie  bij  direct  contact  (met  vlammen)  van  het  muniticartikcl  tn  een 
brand  (op warming  door  si  rating  ),  voorde  120  mm  granaat  met  TNT, 


P  [K/S] 

texpio&re  [s] 

Tuu*,  [XI 

T„end  [X] 

N.v.t. 

81 

947 

382 
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Samenvattend  120  mm 

Uit  eerdere  berekenmgen  [3]  en  laboratoriumexperimenien  met  HMX  houdende 
explosieve  stoffen  is  af  te  lei  den  dat  deze  bij  een  opwarmsnelheid  van  ongeveer 
0,001  °C/s  een  reactie  temperatuur  hebben  bij  182  °C  en  bij  I  °C/s  boven  de  260  °C  met 
een  lineair  logaritmisch  verband,  Deze  waarden  zijn  toegevoegd  in  figuur  4,13a. 

Figuur  4.13a  laat  de  explosietemperaturen  van  RDX,  TNT  en  HMX  als  functie  van  de 
opwarmsnelheid  zien,  RDX  reageert  sneller,  dat  wil  zeggen  het  RDX  loopt  bij  een 
lagere  temperatuur  en  een  kortere  tijd  weg  dan  TNT  of  HMX.  Daarom  gee  ft  RDX  een 
goed  bee  Id  van  de  situatie  in  de  worst  ease. 


a  nn 


T  jwconv  RDX 
- RDX 


C 

o 


■“  200 


150 


100 


0,001 


0,01 

Heating  rate  [°C/s] 


0,1 


0,0001 


Figuur  4.13a  Explosietemperatuur  als  tunc  tic  van  de  opwarmsnelheid  voor  TNT  en  HMX  ten  opzichte  van 


RDX  voor  een  I2t)  mm  granaat. 


Samenvattend  kan  geconeludeerd  warden  dat  voor  de  meest  gangbare  explosieven  in 
conventionele  munide  (RDX*  HMX*  TNT  en  comp  B  ( RDX/TNT)*  Comp  A3  (RDX) 
etc  . )  RDX  de  laagste  waarde  voor  de  explosietemperatuur  geeft  Dit  explosief  gee  ft  bij 
relatief  lage  opwarmsnelheid  van  0*001  °C/s  (3,6  °C/uur!)  pas  een  reactie  boven  de 
150  °C.  Experimented  is  al  gevonden  [3]  dat  een  compositie  waarin  zich  RDX  bevindt 
nooit  veel  hoger  dan  bij  203  QC  zal  reageren  (bij  hoge  opwamisnelheden).  Ook  voor 
veel  moderne  explosieven  waarbij  RDX  en  HMX  als  hoofdbestanddeel  in  een 
polymeren  binder  (HTPB)  zijn  ingebed  (PRXen),  gelden  dezelfde  getallen. 

In  figuur  4. 1 3b  zijn  de  veilige  en  onveilige  gebieden  samengevat  weergegeven.  Ook  een 
nieuw  aangeschafte  155  mm  munitie  van  Rheinmetail  met  de  explosieve  stof  KS  32  vaU 
onder  deze  richtlijnen. 
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250 


Gebi^d  mol  Explosia  gevaar 


Vailig  gebied 


IE  3 


o,oi  o,t 

O  p  warms  nelheid  |*C/s] 


Figuur  4, 1 3b  Samenvaiting  van  de  gebicdcn  waarbij  met  dc  meesi  gangbare  munilie  gevaren  bij  brand 
verwachl  kurrnen  worden  en  waar  met. 


4.2  Berekeningen  aan  kruiten 

4.2  J  Inleiding 

In  de  verschillende  typen  munitie  worden  verschillende  kruiten  gebruikt,  bijvoorbeeld 
single  base,  double  base  en  triple  base  kruiten.  Ook  de  hoeveelheid  en  het  type 
stabilisator  in  bet  kruit  kan  per  artikel  verschillen.  Bij  de  155  mm  granaat  hoort  een 
kardoes  met  kruit.  Figuur  4.14  en  4.15  laien  het  kruit  uit  een  kardoes,  respectievelijk 
een  kardoes  zien. 

Voor  de  120  mm  zijn  twee  typen  mogeSijk.  De  120  mm  mortieren  wrorden  afgevuurd 
met  ringvormig  kruit.  De  120  mm  patronen  worden  door  kruit  in  een  huls  voort- 
gedreven.  Figuur  4.6  en  figuur  4,7  in  de  vorige  paragraaf  lonen  de  120  mm  munitie. 
Deze  voortdnj vende  lading  is  gebaseerd  op  onder  andere  een  double  base  propellant 
met  een  deel  van  de  nitroglycerine  vervangen  door  DEGDN  (diethyleenglycoldimtraat). 

Als  voorbeeld  van  de  35  mm  munitie  wordt  het  brisantgranaatpatroon  DM1 1  met  stalen 
huls  en  schokbuis  DM31 1  in  schakelband  genomen.  De  35  mm  brisantgranaatpatroon 
wordt  gebruikt  in  het  vvapensysteem  van  de  PantserRups  Tegen  Luchtdoelen  (PRTL). 
De  FRTL  is  uitgerust  met  2  computergestuurde  radargeleide  35  mm  snelvuurkanonnen 
en  wordt  gebruikt  voor  de  beseherming  van  onder  meer  bruggen  en  kwetsbare 
installaties  tegen  aanvallen  van  uit  de  lucht.  De  afmetmgen  van  de  35  mm  staan  vermeld 
in  label  4,9,  Figuren  4. 16  tot  en  met  4,18  laten  wat  afbeeldingen  zien  van  de  35  mm 
munitie. 
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Rguur  4. 14  Kruit  till  ccn  white  bag  kardoes  voor  dt  1 55  mm. 


M4C4 


Hguur  4. 1 5  Kardoes  am dc  granaat  van  do  155  mm  weg  le  schieteiL 


Tabe  14.9  A fme t  i  rtge  n  3  5  mm  pal  rk i  n . 


Af  meting 

in  mm 

lengte  patroon 

387 

lengte  granaat 

188 

lengte  schokbuis 

75 

diameter  granaat 

35 

diameter  hulsbodem 

54(6 
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Figuur  4. 16  Links:  b  ri  sa  nr  g  ranaai  pam>on  DM  1 1 ,  rectos:  ops  lag  van  de  patronen  op  pallet 


VS  9-447 
6  07.01 
WrtdS 


1  fctxMtxMs  QWP  t  (n  tfb  31 

2  grMnnt&ch**m 
S  xpnngifftJtnQ 

4  gtfimttwxl 

5  1XK*n**H; 

0  iadUnp 

l  putwonhute 

6  QntsfkfaQgtioG  f/i*  afb  4) 


Figuur  4. 17  Links:  aanzicht  van brisantgranaatpairoon  DM  II,  rectos:  dtvorsnede  van  brisantgranaatpairoon* 
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Figuur  4. 1 S  Patronen  in  een  schakelband 

In  paragraaf  3.2  werd  a  I  toegelieht  dat  de  invoerpara  meters  voor  kmit  iastig  te  bepalen 
zijn.  Om  deze  redenen  is  ervoor  gekozen  om  voor  de  voorgestelde  7  koppels  van 
invoerparameters  (kruit  best  case  1,  kruit  best  case  2,  kruit  worst  case  1 .  kmit  worst 
case  2,  KB6954,  kruit  C  en  kruit  M)  de  tijd  en  waar  mogelijk  de  temperatuur  te  bepalen 
voor  een  cilindrische  vorm  met  een  diameter  van  120  mm.  Tijdens  de  l-dimensionale 
berekeningen  bleek  dat  de  bandbreedte  van  het  kruit  (best  en  worst  case)  zoals 
voorgesteld  niet  de  explosietijden  en  temperaturen  van  de  voorbeeldkruiten  ( KB6954. 
kruit  C  en  M)  bevatte.  Blijkbaar  zijn  de  voorgestelde  parameters  niet  breed  genoeg  om 
deze  bandbreedte  op  te  siellen.  Ter  vergelijking  naast  de  resultaten  uit  het 
l-dimensionale  model  is  de  omzetting  van  de  stabilisator  berekend,  Deze  omzeUing  is 
temperatuursafhankelijk  en  derhalve  is  bij  verschillende  opwarmsnelheden  deze 
converste  bepaald. 

Om  de  resultaten  in  label  vorm  aan  te  leveren  zal  te  onoverzichtelijk  worden. 

Daarom  zijn  de  resultaten  in  een  grote  grafiek  weergegeven  (figuur  4.19)  zodat  een 
bandbreedte  zichtbaar  wordt  waarbinnen  de  kmiten  reageren  bij  verschillende 
opwamisnelheden*  Figuur  4. 19  be  vat  de  op  warming  via  direct  contact  van  bet  kruii 
(groene  curves)  en  de  temperaturen  waarbij  de  stabilisator  verbruikt  is  (rode  curves), 
Om  de  bandbreedte  te  valideren  zijn  er  experimenten  uitgevoerd.  in  de  voigende 
paragraaf  worden  de  resultaten  van  de  experimenten  met  het  KRAI 386  kruit  gegeven. 
De  resultaten  zijn  al  opgenomen  in  onderstaande  figuur. 
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Qpwarmsnelheid  [K/s] 


—  Temp@ra.luur  direct  contact  -  KB6954 
Temps raluur  direct  contact  -  kruil  C 

’*  Temperaluur  direct  contact  -  kruil  M 
-»  Temperatuur  direct  contact  -  best  esse  1 

—  Temperatuur  direct  contact  ■  best  case  2 
— Temperatuur  direct  contact  ■  wore!  case  t 

TempeTatuur  direct  contact  -  worst  case  2 
Temperatuur  omgeving  kruil  KB6954 
-m-  Temperatuur  omgeving  kruit  C 
Tflm  peraluur  omgevrng  krurt  M 

—  Temperatuur  omgevmg  kruit  best  case  1 
Temperatuuromgeving  kruit  best  case  2 
Temperatuur  orngeviirg  kruit  worst  case  1 
Temperatuur  omgeving  knjit  worst  case  2 

"r^Enpenment 


Figuur  4. 1 9  Vergelijking  van  temperatuur  waarbij  wegloop  plaat.s  vindt  hi  j  opwarming  door  direct  contact 
(groen)  en  consumplie  van  tie  stabilisator  (rood). 


4. 2. 2  Expe  rimen  ten  m  et  kruit  va  n  G  ait  ling  30  m  m  Gun  ( K  RA  1386) 

Wereldwijd  zijn  er  niet  zoveel  opwarmexperimenten  uitgevoerd  met  kruiten  om  te 
bepalen  bij  welke  temperatuur  deze  reageren.  Om  een  indicatie  te  hebben  van  de 
responsie  van  kruiten  is  er  nog  een  testserie  uitgevoerd  met  de  SCB  [14]  test  met  het 
Nitrocellulose  kruit  (single  base)  van  Goalkeepermunltie  (30  mm).  Figuur  4.20  laat  een 
SCB  vaatje  zien  gevuld  met  dit  kruit. 


Figuur  4.2U  N itrocel lu losekru i t  in  een  SCB  vaatje  (bovenaanzicht) 


In  figuur  4.21  is  een  typisch  resultaaf  te  zien  na  het  experiment.  Hr  heeft  duidelijk  een 
explosie  plaats  gevonden.  In  tabel  4. 10  is  het  meetoverzicht  gegeven  van  deze  testserie. 
De  resultaten  zijn  eveneens  in  het  overzicht  van  figuur  4.19  terug  te  vindem 
De  experimenteie  gegevens  vallen  binnen  de  range  van  alle  berekende  resultaten. 
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label  4  10  Resultatcn  van  de  SCR  lestserie  mei  KRAI 386  kruil  van  do  Goalkeepermunitie 


Opwarmsnelheid  [KJs] 

Temperatuur  van  de  wand  [°C] 

0.06 

166.  169 

0,01 

_ 152 

0,0033 

142 

piguur  4.21  Typische  responsie  van  KRA  13K6  kruil  in  een  opwarmtesl. 


Samenvallend  kan  het  volgende  geeoncludeerd  worden,  Uil  de  here  ken  ingen,  maar  ook 
uit  de  experimentele  gegevens  blijkt  dat  kruiten  eerder  reageren  bij  opwarming  dan  de 
explosieven  die  in  veel  gangbare  munitie  zit.  De  geextrapoleerde  waarde  voor  een 
opwarmsnelheid  van  0.001  °C/s  van  het  kruit  in  de  experimenten  is  ongeveer  125  °C  en 
de  berekening  met  de  worst  case  is  dil  bij  een  opwarmsnelheid  van  0,00loC/s  zelfs 
105  °C.  In  hoeverre  deze  waarde  echt  reeel  is.  is  niei  preeies  in  te  schatten  maar  deze 
waarde  kan  zeker ook  niet  uitgesloten  worden.  Hierbij  is  het  scheidingsgebied  voor 
gevaar  en  geen  gevaar  een  rechte  lij n. 

200 


o  1M- 


50- 


Gebied  met  Explore  gevaar 


Veilig  gebied 


- ' — ' — * . .  * — ' — . . . - — ' — . 

1E3  CGI  0J  t 

Opwarmsnelheid  [°C/s] 


Hguur  4,22  Samenvattmg  van  dc  gebieden  waarbij  met  de  tneesl  gangbare  krniten  gevaren  bij  brand 
verwachl  kunnen  worden  en  waar  niet. 
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5  Afhankelijkheid  van  warmtegeleidbaarheid,  warmte- 

capacitcit  en  dichtheid  op  de  explosictemperatuur  en  -tijd 


Om  de  afhankelijkheid  van  de  thermische  geleidbaarheid  h  de  wafnttecapaciteit  en  de 
dichtheid  van  het  TNT  op  de  cook-off  tijd  te  onderzoeken,  zijn  met  de  overige 
eigensehappen  uit  [1],  een  opwarmsnelheid  van  0 A  K/s  en  een  diameter  120  mm  de 
cook-off  tijden  berekend  een  range  van  waarden.  Deze  waarden  zijn  de  litemtuur- 
waarden  van  deze  parameters,  een  50%  kleinere.  een  25%  kleinere  waarde,  een  25% 
grotere  en  een  50%  grotere  waarde,  Tabel  5.1  tot  en  met  5.3  geven  overzichten  aan  van 
de  gebruikte  w  aarden  voor  X.  p  en  Cp  en  de  berekende  cook-off  tijden. 


Ta bo  15.1  Cot  ik-off  tijde n  bij  ve  rsehi  1  le  ndo  waarde  n  voor  de  ge  lo  idbaa  rhe  id . 


A 

Afwijking  in  A 

texplcxte 

Afwijking  in 

[W/m-K] 

[%i 

[s] 

[%] 

0,1046 

-50 

4658 

-0 

0,1569 

-25 

4888 

-4 

0,2092 

0 

5080 

0 

0,2615 

+25 

5241 

3 

0,3138 

+50 

5377 

6 

Overige  parameters  uii  Elf  opwarmsnelheid  is  0,1  K/s, 

Tahel  5.2 

Cook-off  tijden  bij  verschillende  waarden  voor  de  dichtheid, 

P  [kg/m3] 

Afwijking  in  p  f%] 

t&xplsie  [Sj 

Afwijking  in  tBlptoK„  [%] 

785 

-50 

4237 

-17 

1177.5 

-25 

4704 

-7 

1570 

0 

5080 

0 

1962.5 

+25 

5402 

6 

|  2355 

+50 

5687 

12 

Overige  parameters  uii  |  i  f  opwarmsnelheid  is  O.  I  K/s 


Tabel  5.3  Cook-off  tijden  bij  verscbillende  waarden  voor  de  warmlecapacitcil. 


Overige  parameters  uii  1 1  f  opwarmsnelheid  is  0,1  K/s, 
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2 

QJ  ^ 

CL 

~  n> 

O  8 

°  I 

trj 


5 


(  jguur  5 . 1  Sproiding  van  dc  cook-off  Ujden  bij  afwijkende  warmtcgeleidbaarhdd.  dichlhcid  on 
warmtccapac  iteit . 


Wanneer  de  variaties  in  de  geleidbaarheid,  de  dichtheid  en  de  warmtecapaciteit  in 
figuur  5. 1  bekeken  worden,  blijkt  dal  de  variatie  in  cook-off  tijd  het  minst  is  bij  een 
variatie  in  de  geleidbaarheid  (X).  De  variatie  in  warmtecapaciteit  leidt  tot  de  grootste 
variatie  in  cook-off  tijd. 


Bovenstaande  resultaten  laten  zien  dat  de  resultaten  van  het  mode)  niet  sterk  afhankelijk 
zijn  van  de  parameters  X,  p  en  Cp.  In  praktijk  kunnen  deze  parameters  goed  binnen  de 
bovenstaande  nauwkeurigheid  bepaald  worden.  De  resultaten  voor  de  responsie  zullen 
dan  ook  betrouwbaar  zijn. 
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6  Vergelijking  tussen  explosieven  met  hoog  en  laag  RDX 
gehalte 


Het  vorige  hoofdstuk  behundelde  een  gevoeligheidsstudie  van  de  parameters  in  het 
model  In  dit  hoofdstuk  is  voor  een  case  de  uitkomst  van  het  model  vergeleken  met 
waarden  uit  de  literatuur. 

Scholtes  et  ak  [9]  hebben  de  cook-off  temperatuur  (de  temperatuur  van  de  wand  op 
moment  van  exploderen  van  het  monster),  van  een  monster  PRXN-109  (een  veel 
gebruikt  explosief  voor  gevechtsiadingen)  experimented  bepaald.  PBXN-109  bestaal 
uit  64  wt%  RDX,  16  wt%  binder,  20  wt%  aluminium.  In  het  model  is  dit  gesimuleerd 
als  een  massieve  eilinder  met  dezelfde  diameter  met  de  parameters  van  RDX  met  een 
gecorrigeerde  waarde  voor  de  reactiewarmte,  namelijk  64  wt%  in  plaats  van  85  wl%, 
RU8I  hestaat  voor  85  wt%  uit  RDX  met  15  wt%  binder 


Tabs  16.1  Sto  fa fhanke  li  j  ke  e i  go  nschappe n  van  de  e  xp  los  ie ve  n  gcb  ru  iki  i  n  de  he  re  ke  ni  nge  n . 


Parameter 

Symbool 

Eenheid 

RU81 

PBXN-109 

Warmtegeieidbaarheid 

A 

W/(m-K) 

0.25 

0.25 

Dichtherd 

P 

kg/m3 

1580 

1580 

Warmtecapaciteit 

Cp 

J/(kg  K) 

1360 

1360 

Reaetiewarmte 

Q 

J/kg 

20888900.85' 

2088890  0. 64 r 

Pre-exponentiele  coefficient 

Z 

1/s 

3,16228-IQ18 

3,16228-1 01fl 

Activeringsenergie 

Ea 

J/mol 

199000 

199000 

Explosietemperatuur 

Texplceia 

*C 

203 

203 

{  De  compos i tie  hestaat  voor  deze  fraetie  uit  het  explosief 


Tabd  6 . 2  Re  ke  npar ameiers  voor  het  mode  l 


Parameter 

Symbool 

Eenheid 

Waarde 

Opwarmsnelheid 

P 

K/s 

0,08333 

Warmteoverdrachtscoefficient 

n 

W/(m 2K) 

8 

Mogelijkheid  tot  opwarming 

optie  1 

- 

1 

Op  warming  via  vrije  convectie  of 
direct  contact 

optie  2 

* 

1 

Mogetijkheid  voor  point  of  no  return 

optie  3 

- 

0 

Mogeiijkheid  voor  de  stralingswarmte 

optie  4 

* 

0 

Straai  van  mogeltjke  bore,f 

Rb 

m 

0 

Straal  van  explosief' 

RP 

m 

0,0233 

Straai  van  liner1 

Ri 

m 

0,031725 

Tijdstap 

t 

s 

varieert 

Starttemperatuur 

Tstari 

K 

293.15 

O  mg  ev  i  n  gste  mpe  rat  u  u  r 

"fcing 

K 

293,15 

J  In  de  berekeningen  is  aangenomen  dat  de  cilinder  uit  alleen  RDX  besiaat  cn  geen  'bore'  of  liner  beval. 


Voor  het  RDX  geldt  dat  zodra  de  temperatuur  in  het  explosief  boven  de  203  °C  komt  de 
wandtemperatuur  genome  n  wordt  als  temperatuur  waarbij  de  cook-off  plaats  vindt. 

In  tabel  6.3  staan  de  resultaten  van  de  experimenten  uit  refencie.  [9]  en  van  de 
berekeningen.  De  resultaten  zijn  berekend  met  de  opwarming  van  de  wand  by  direct 
contact  aangezien  btj  het  experiment  een  verwarmings mantel  de  wand  direct 
verwarmde. 
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label  6.3  Lineaire  opwarming  RDX-PBX. 


Tijd  [S] 

1"wand 

Experimenteel  RDX  PBXN-109 

18953 

179 

Experimenteel  i-RDX  PBXN-109 

19764 

182 

Modal  64%  RDX 

19415 

182 

Model  85%  RDX 

19135 

179 

De  berekende  waarden  van  de  tijd  tot  explosie  en  de  temporal uur  komen  overeen  met  de 
gemeten  waarden.  Bovenstaande  residtaten  laten  zien  dat  een  lagere  concentratie  RDX 
in  combi natie  met  aluminium  ook  goed  benaderd  wordt  met  het  model.  Hei  model  kan 
dan  ook  gebruikt  worden  voor  explosieven  die  minder  RDX  bevatten  en  meer  niet- 
energetisch  maternal. 
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7  Bulkopslag  kruit 

De  uitgevoerde  berekeningen  geven  de  temperatuur  aan  waarboven  bet  an  ike  l  reageert 
bij  op  warming  van  de  afzonderlijke  art  ike  ten.  De  invloed  van  bulkopslag  op  de 
reaetietijd  en  temperatuur  is  in  dit  hoofdstuk  bekeken,  Een  volledige  behandeling  van 
deze  theorie  is  gegeven  in  bijlage  A. 

7.1  Veiligc  opslag  temperatuur  voor  kruit 

De  kritische  omgevingstemperatuur  voor  een  stapeling  van  kruit  is  bepaald  met  hulp 
van  een  simulatiemodel  door  Lont  [1  ].  Een  stapeling  van  twee  lagen  kruit  in  de  hoogte 
van  5  bij  5  vaten  is  gesimuleerd.  De  warmtestroming  in  de  ruimte  tussen  de  vaten  is 
gemodelleerd  met  behulp  van  een  effective  warmtegeleidingcoefficient  van 
1  W/(m-K).  Het  kruit  zit  opges  lagen  in  een  fiberdrum  met  de  volgende  a  f  met  ingen: 
hoogte  =  630  mm,  diameter  =  355  mm  met  wanddikte  =  2,3  mm. 


Tabe  17.1  Kritische  te mperat u re n  van 

1  vai  kmii  cn  een  stapeling  van  vaten  kruit. 

KB6954 

Wh  1  drum*.  T1  [°C) 

82,0 

T*,.,**  2*5*5  drums*.  Ts  [  C] 

68,6 

Verschil  (Tl-Ts)  [°C] 

13,4 

*  Drum:  diameter  =  355  mm,  heigih  =  630  mm. 


De  simulatie  laat  zien  dat  de  kritische  temperatuur  van  bulkopslag  lager  ligt  dan  de 
kritische  temperatuur  van  de  afzonderlijke  drums  kruit.  zie  label  7.1 .  De  stapel  had  een 
10  °C  tot  15  °C  lagere  kritische  omgevingstemperatuur  dan  een  enkel  vat.  De  laagste 
kritische  omgevingstemperatuur  uit  de  simulates  is  circa  50  °C 

7*2  Ttr  volgens  Boddington  simulatie  kubus 

Een  cweede  bulkopslag  is  gesimuleerd  met  een  model  beschreven  in  bijlage  A.  De  bulk¬ 
opslag  is  gesimuleerd  als  een  kubus  van  0,5  x  0.5  x  0.5  m,  Met  de  gegevens  van  het 
kruit  met  naam  KB6954  wordt  met  de  zogenaamde  formule  van  Kamanetski  [3]  een 
waarde  van  de  kritische  temperatuur  van  78.0  °C  gevonden  terwijl  met  een  eindige 
elementen  methode  ABAQUS  een  waarde  van  77.9  °C  wordt  gevonden.  Om  het 
temperatuur  profiel  te  berekenen  wordt  de  zogenaamde  methode  van  Boddington 
gebruikt,  Samenvattend  zegt  deze: 

De  methode  van  Boddington  is  geldig  in  het  Frank-Kamenetskii  regiem.  Uit  de 
berekeningen  blijkt  dat  een  5C  orde  benadering  nodig  is.  Bij  een  T  orde  benadering 
wijken  de  kritische  temperaturen  teveel  af  in  verge lij king  met  de  berekeningen  uit 
ABAQUS.  Een  5C  orde  benadering  gee  ft  betere  resultaten.  In  label  7.2  staan  de  kritische 
temperaturen  weergegeven.  De  afstand  van  het  nulpunt  tot  het  punt  van  de  berekening 
is  weergegeven  als  p,  De  gebruikte  parameters  zijn  opgenomen  in  tabel  7.3. 

De  berekeningen  zijn  uitgevoerd  voor  een  kubus  met  een  ribbe  van  0,5  m. 
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'l  abel  12  Temperatmuprofiel  kuhus  vtx>r  5e  orde  benadering  van  Boddingtnn. 


p 

T  (K) 

T  (°C) 

P 

T  (K) 

T(’C) 

0,00 

363,02 

89,87 

0,52 

350,38 

85,23 

0,04 

362,99 

89,84 

0,56 

357.78 

84,63 

0,08 

362,89 

89,74 

0,60 

357,17 

84,02 

0,12 

362,72 

89,57 

0,64 

356,55 

83,40 

0,16 

362,50 

89,35 

0t68 

355*93 

82,78 

0,20 

362,21 

89,06 

0,72 

355,30 

82,15 

0,24 

361,87 

88t72 

0,76 

354,67 

61,52 

0,28 

361,48 

88,33 

0,80 

354,04 

80,89 

0,32 

361,05 

87,90 

0,84 

353,40 

80,25 

0,36 

360,57 

87,42 

0,88 

352,77 

79,62 

0,40 

360,06 

86,91 

0,92 

352,13 

78,98 

0,44 

359,52 

86,37 

0,96 

351,48 

78,33 

0,48 

358,96 

85,81 

1,00 

350,82 

77,67 

Tabol  7.3 

Parameters  gebruiki  in  de  herekenmg  volgens  de  methods  van 

Boddingti>n, 

K 

J/(msK) 

0,08 

P 

kg/m3 

840 

Q 

J/kg 

1 

Z 

1/s 

5.782*1 0Te 

Ea 

J/mol 

1,227  "10s 

Bovenstaande  waarden  zijn  de  fysische  eigenschappen  van  het  kruit  KB  6954.  zoals  ook 
gesimuleerd  als  kruit  B  in  [JO]. 


Duidelijk  is  voor  hulkopslag  dat  de  kritische  temperatuur.  waarbij  een  kruit  langzaam 
begint  te  reageren.  a  I  rond  de  70  graden  kan  (out  of  area  situaties)  liggen.  afhankelijk 
van  het  kruit,  Als  deze  opslagtemperatuur  lang  gehandhaafd  blijft  kan  dit  leiden  tot  een 
wegloopreactie  van  het  kruit  en  dus  een  gevaarlijke  situatie.  Theoretisch  duurt  het  zo'n 
230  dagen  voordat  dit  echt  wegloopt,  echter  de  kwaliteit  van  het  kruit  zal  hierdoor 
drastisch  achteruitgaan,  Bij  hogere  temperaturen  zal  dit  nog  sneller  gaan. 
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8  Conclusies  en  aanbevelingen 


De  grootste  oorzaak  bij  ongelukken  met  munitie  is  brand.  Met  name  conventionele 
muni  tie  geeft  in  een  brand  of  in  de  nabijheid  van  brand  grote  gevaren  voor  omgeving  en 
personeel  Hoe  te  handelen  bij  brand  wordt  ondermeer  omschreven  in  de  richtlijnen  van 
MP  40-21  en  de  brandweeropleiding  Brandmeester,  In  deze  studie  is  getracht  enkele 
methoden  aan  te  geven  om  deze  richtlijnen  aan  te  passen.  Na  enkele  bezoeken  van 
munitie  opslag  bij  Defensie,  een  literatuuronderzoek  naar  de  chemische  componenten 
van  belang  bij  munitie  in  een  brand  zijn  de  voorhanden  rekenspreadsheet  en  een 
I-dimensionale  code  aangepast  om  berekeningen  uit  te  voeren  aan  35,  1 20  en  155  mm 
munitie  die  voorradig  is  bij  de  KL.  Behoudens  explosteve  stoffen  is  er  ook  gekeken 
naar  voortdrijvende  ladingen  van  dit  soort  muni tiety pen.  De  volgende  conclusies 
kunnen  getrokken  worden. 

Het  wordt  aanbevolen  om  wel  camera's  en  temperatuursensoren  te  plaatsen  zodat  op 
afstand  mogelijke  opwarming  geidentificeerd  kan  worden  en  de  aard  en  loeatie  van 
deze  opwarming  bepaald  kan  worden.  Na  aflezing  van  de  geregistreerde  temperatuur  en 
de  daaruit  afgeleide  opw  armsnelheid  kan  eventueel  besloten  worden  of  er  nog  actie 
genomen  wordt  om  de  brand  te  hlussen  of  dat  er  op  een  afstand  van  1000  m  gewacht 
wordt  totdat  een  eventuele  explosie  heeft  plaatsgevonden. 

Omeen  benadering  te  krijgen  voor  de  tijd  die  men  nog  heeft  om  in  te  grijpen  wanneer 
munitieartikelen  met  explosieven  opwarmen  kan  de  worst  case  grens  van  RDX 
aangehouden  worden.  Voor  gevechtsladingen  ligl  deze  grens  op  de  lijn  van  150  °C  bij 
een  opwarming  van  0,001  K/s  lineair  naar  200  °C  bij  een  opwarming  van  0,1  K/s  en 
opwarmsnelheden  daar  boven.  De  tijd  die  daarmee  gepaard  gaat  is  afliankelijk  van  de 
temperatuur  op  moment  van  aflezen  en  de  opwarmsnelheid,  Kruiten  reageren  bij  lagere 
temperaturen.  Voor  artikelen  met  kruit  ligt  deze  grens  lager,  van  105  °C  naar  150  °C  bij 
snelheden  van  0.001  K/s  lot  l  K/s,  Voor  munitie  in  een  directe  brand  is  uiterste 
voorzichtigheid  geboden.  Tijden  van  minder  dan  1  minuut  zijn  met  vreemd. 

Om  een  beter  overzicht  van  de  kruiten  te  krijgen  zal  een  grotere  variatie  van  ‘snel 
reagerende  kruiten'  bekeken  moeten  worden. 

De  meest  gangbare  explosie ven  in  conventionele  munitie  (RDX.  HMX,  TNT  en  comp 
B  (RDX/TNT),  Comp  A3  {RDX  et  cetera),  kunnen  benaderd  worden  met  de  berekening 
van  RDX.  Dit  geeft  namelijk  de  laagste  waarde,  Dit  geldt  ook  voor  veel  moderne 
PBXen  met  als  hoofdbestanddeet  RDX,  HMX  of  TNT  en  bindersystemen  als  HTPB. 
Ook  bijvoorbeeld  de  moderne  155  munitie  van  Rheinmetall  (met  KS32)  valt  binnen 
deze  normen,  net  als  munitie  met  het  veelgebruikte  'blast'  explosief  PBX-N109, 

Het  ontw'ikkelde  model  om  de  temperaturen  te  bepalen  is  niet  sterk  afhankelijk  van  de 
parameters  Xf  p  en  Cp,  In  praktijk  kunnen  deze  parameters  goed  binnen  de 
bovenstaande  nauwkeurigheid  bepaald  worden,  De  resultaten  voor  de  responsie  zullen 
dan  ook  betrouwbaar  zijn. 
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Voor  bulkopslag  ligt  de  kritische  temperatuur  waarbij  een  kruit  langzaam  begint  te 
reageren  al  rond  de  70  graden,  afhankelijk  van  bet  kruit.  Als  deze  opslagtemperaiuur 
lang  gehandhaafd  blijft  kan  dit  leiden  tot  een  wegloopreactie  van  bet  kruit  en  dus  een 
gevaarlijke  situatie,  Tbeoretisch  duurt  bet  zo' n  230  dagen  voordat  dit  ecbt  wegloopt, 
Fxhter  de  kwaliteit  van  bet  kruit  zal  hierdoor  drastisch  achteniitgaan.  Bij  hogere 
temperaturen  versnelt  dit  proces  en  zal  de  levensduur  snel  aehteruit  gaan.  Hierdoor 
kunnen  eveneens  gevaarlijke  situaties  ontstaan. 
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A  Berekening  thermische  explosickubus 


A.l  Inleiding 


De  theorie  van  een  thermische  explosie  is  in  de  eerste  he  I  ft  van  vorige  eeuw  ontwikkeld 
door  Semenov'  en  Frank- Kamenetskii2.  Het  verschil  tussen  beide  benaderingen  is  dat 
bij  Semenov  de  inwendige  temperatuur  van  het  lichaam  constant  is  en  een 
temperatuurs  verschil  bestaat  tussen  lichaam  en  omgeving,  terwijl  bij  Frank- 
Kamenetskii  er  een  temperatuursgradient  bestaat  in  het  lichaam  terwijl  de  oppervlakte- 
temperatuur  gelijk  wordt  gesteld  aan  de  omgevings temperatuur.  De  Frank- Kamenetskii 
benadering  is  het  meest  realistisch  in  het  geva!  van  stoffen  met  een  lage  warmte- 
geleidingcoefficient.  zoals  explosieve  stoffen  en  de  meest  recente  theorieen  zijn  daarom 
op  Jit  principe  gebaseerd.  De  meeste  theorieen  zijn  gebaseerd  op  de  Frank- Kamenetskii 
parameter  5  die  gedefinieerd  is  voor  de  drie  geometrisch  simpele  configuraties  waarin 
de  temperatuur  een  functie  is  van  slechts  ddn  coordinaat:  de  bol,  de  oneindige  cylinder 
en  de  oneindige  vlakke  plaat.  Deze  parameter  wordt  gegeven  door: 


oQAEa2e~E,Kr* 

ART'; 


(A.l) 


waarin  a  de  straal  of  halfbreedte  is.  TA  de  omgevingstemperatuur,  R  de  gasconstante, 
E  de  activeringsenergie.  o  de  dichtheid,  Q  de  reactiewarmte  per  massaeenheid,  A  de 
Arrhenius  factor  en  X  de  warmtegeleidingscoefficient. 


Voor  andere  geometrieen  is  het  in  principe  niet  mogelijk  oni  de  temperatuurprofielen 
als  functie  van  slechts  e£n  coordinaat  te  geven.  Het  is  in  dat  geva]  ook  met  mogelijk  een 
algemene  analytische  oplossing  te  geven.  Numerieke  oplossingen  voor  gedefinieerde 
gevallen  zijn  uiteraard  wel  mogelijk.  maar  het  bepalen  van  de  kritische  condities  is  een 
vrij  omslachtige  operatie  en  levert  weinig  inzicht  in  de  thermische  eigensc happen  op, 
Om  hieraan  het  hoofd  te  hieden  is  door  Boddington  et  al/  een  benaderende  algemene 
oplossing  gegeven  voor  een  bepaalde  klasse  van  voorwerpen  voor  het  Frank- 
Kamenetskii  regiem,  Deze  benadering  is  van  toepassing  op  voorwerpen  die  een 
symmetriecentrum  hebben  en  waarbij  een  lijn  vanuit  het  centrum  het  oppervlak  slechts 
ednmaa!  snijdt.  Hierbij  wordt  de  oplossing  gegeven  in  termen  van  een  enkele  variabele 
p  die  is  gedefinieerd  als  voIgt:  p(P)  is  de  afstand  van  het  punt  P  in  het  lichaam  tot  het 
symmetriecentrum  van  dat  lichaam,  gedeeld  door  de  afstand  van  het  centrum  O  tot  het 
oppervlak  S  in  de  rich  ting  OP.  Op  deze  manier  wordt  voor  het  lichaam  een  'equivalente 
bol’  gecreeerd  waarop  de  thermische  explosie  theorie  wordt  toegepast. 
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Algemene  oplossing  volgens  Boddington 


Bij  gebruik  van  de  volgende  definities: 


0=  E1{T- 
RT‘ 

ta) 

(A.2) 

met  0O  —  0  in  centrum. 

p  =  dimensieloze  afstand  (0  in  centrum;  1  op  grensvlak) 

1  1  redd) 

7f=  7^ 

(A3) 

(R<>  is  de  afstand  r  van  het  centrum  tot  het  grensvlak,  gemiddeld  over  4k  steradialenk 
wordt  de  algemene  oplossing  gegeven  door: 

0  =  0O  -  2  In 

^<1=0  J 

(A. 4) 

«tk'2 

(A. 5) 

S(R„)=\2U  +  l)  y2 

K  0/  3  X? 

(A.6) 

Hierbij  is: 

£(*o)=J*o 

(A. 7) 

oQAEe-E,RT ' 

met:  Y  ~  ■  - 

MiT: 

(A.8) 

Deze  definitie  sluit  aan  bij  de  definitie  van  die  van  Frank- Kamenetskii  gegeven  in  ( A  J  ) 

2 

waarbij  geldt:  8|-k  —  8(a)  =  y  dC  -  — -<?(/?„  )  met  a  =  halfbreedte,  straal.  et  cetera. 

R'o 

Door  Boddington  wordt  de  dimensieloze  temperatuur  0  dus  uitgedrukt  a  Is  een 
machtreeks  in  yp'  waarbij  de  coefficienten  b  worden  gegeven  door: 

K=b,=\ 

(n  +  lX2n  +  j  + 1)6,*,  =^{n-s){As  ~2n  +  3- j)bjnbn_s  {n  >  l) 

.f= 0 

D  2 

met:  j  =  3-^-1  (A.  10) 

a' 


j  is  de  zogenaamde  vormparameter  die  bepaald  wordt  door  de  vorm  van  het  voorwerp, 
terwijl  de  parameter  y  afhankelijk  is  van  0„  en  j  en  bepaald  wordt  door  de  rand- 
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voorwaarden.  Volgens  Boddington  geldt  dat  als  de  reeks  afgekapt  wordt  bij  n=5  de  tout 
ten  opzichte  van  de  exacte  oplossing  zeer  klein  is,  terwijl  ook  bij  afkappen  van  de  reeks 
bij  n=2  de  benadering  nog  zeer  good  is. 


A. 2 


Vergelijking  met  analytische  oplossing 


Door  onder  andere  Bowes-*  is  een  analytische  afleiding  gegeven  van  het  temperatuur- 
profiel  voor  het  Frank*  Kamenetskii  regiem  voor  de  gevallen  waarbij  dit  exact  oplosbaar 
is  (oneindige  vlakke  plaat,  oneindige  cylinder,  bol).  Voor  de  oneindige  vlakke  plaat  is 
dan  de  oplossing: 


( 


0  =  0,  -  2  In  cosh 


Se 9" 


\ 


V  2 


(A.  11) 


©rt/2 

_  € 


{AM) 


of: 


0  =  0„ 


2  In  cosh 

v 


2 


(A.  13) 


Om  nil  te  bekijken  hoe  goed  Boddington  en  Bowes  overeenkomen  in  dit  geval  is  de 
aigemene  oplossing  van  Bowes  gesehreven  als  een  maehtreeks  in  yp"  en  vergelijken 
met  de  reeks  van  Boddington.  Voor  dit  geval  geldt  dat  j  =  0  dus  voIgt  uit  (A.6)  dat 
6  =  4y/Xj2  Wanneer  dil  wordt  gesubstitueerd  in  vergelijking  (Ad 3)  van  Bowes  en  de 
cosh  als  een  maehtreeks  wordt  gesehreven  is  het  resultant: 

0  =  ©0  -  2  In  [  1  +  yp3  +  1/6  (yp2)'  +  1/90  (yp2/  +  0.000362  (yp2)4  +  ...  (A.  1 4) 

Uit  de  recur  sieve  relaties  (A.9)  voor  de  coefficienten  bj  van  Boddington  volgen  in  dit 
geval  voor  de  waarden  van  b,:  b0  =  1 :  b|  =  l ;  bs  &  1/6;  b*  -  1/90:  b4  =  0.000397, 
waardoor  de  oplossing  van  Boddington  er  als  volgt  uitziet: 

0  =  0„  -  2  In  [1  +  yp2  +  1/6  (ypY  +  1/90  (yp2/  +  0.000397  (yp2/  +  (A.  15) 

De  twee  oplossingen  blijken  dus  in  dit  geval  overeen  te  komen  voor  een  ftink  aantal 
termen  in  de  maehtreeks.  pas  in  de  orde  term  treedt  een  klein  verse  hi  1  op,  zodat  het  er 
op  lijkt  dat  de  theorie  van  Boddington  inderdaad  de  juiste  resultaten  oplevert  voor  de 
exact  analytisch  te  berekenen  gevallen,  Dit  geeft  vertrouwen  voor  de  toepasbaarheid 
van  deze  theorie  in  andere,  meer  eomplexe  gevallen. 


Kritische  waarde 


Uit  de  thermische  explosie  theorie  volgt  dat  de  kritisehe  toestand  w'ordt  bereikt  als 
geldt: 


dS 

o 


=  0 


(A.  16) 


=  1 


ofwel: 


(A.  17) 
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In  geval  we  ons  beperken  tot  een  2^  orde  benadering.  die  volgens  Boddington  a I  goede 
resultaten  meet  opleveren,  resulteert  dit  in  de  volgende  uitdrukkingen  voor  de  kritische 
w  a  a  rde  n  i  v  a  n  de  d  i  v  erse  groo  thede  n : 

y=(Arl -ij[i-*2(x, -i)] 

(A.I8) 

X,,r=7  0  +  7) 

4 

(A. 19) 

0m  =  2  Inf  2  0  +  7)1 

(A.  20) 

^,=^0  +  3X;+7) 

(A. 2! ) 

*„(*>- *Vt6) 

(A.  22) 

A.3  Toepassing  op  kubus 

De  vormpara meter  j  kan  direkt  via  vergelijkingen  (A.3)  en  (A.  10)  berekend  worden. 
Voor  een  groot  aantal  geometrieen  is  in  appendix  2  van  Boddington  deze  berekening  al 
uitgevoerd.  Voor  de  kubus  geldt: 

=  1  +  —  ^3  Ro=l.  19447a  (a  =  halve  as  van  kubus)  (A.23) 

R;  JC 

Hieruit  volgt  dat  j  =  3.2803  en  zijn  de  coefficienten  b;: 

bn  =  bi  =  I;  bi  =  -0.18154;  b,  =  0,0646 12;  b4  =  -0.030556;  bs  =  0,017192 

A  Is  we  de  2C  orde  benadering  gebruiken  worden  de  kritische  waarden  voor  de  kubus 
gegeven  door: 

XLer  =  2,570 
0,^  =  1.8878 

ycr  =2,0175 

$M  =  3,6654  (=>  Sja)=  2,5690) 

Dit  pas  sen  we  nu  toe  op  een  kubus,  met  een  ribbe  van  0.5  m  la  —  0,25  m}  en  we 
gebruiken  de  volgende  waarden  voor  de  parameters  van  de  explosieve  stof: 

E  =  1.227.10s  J/mol;  A  =  5.782.1016  s  ';  Q=  I  J/kg;  a  =  840  kg/m3;  A  =  0,08  J/s/nVK; 
R  =8.314  J/mol/K 

Hieruit  volgt  voor  de  kritische  waarden: 

Ro  =  0.29862  m 

TA>cr  =351,25  K  =  78.10  °C 

To.cr  =  367.03  K  =  93.99  °C 
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Het  temperatuurprofiel  kan  nu  ook  berekend  worden  en  de  resultaten  zijn  gegeven  in 
tabel  A.  1  en  figuiir  A.  1 .  Ter  vergelijking  is  tevens  in  figuur  A.  1  de  uitkomst  van  een 
numerieke  simulatie  getekend,  uitgevoerd  met  ABAQUS  voor  dezelfde  materiaal- 
parameters  en  met  gebruik making  van  een  rooster  van  9  x  9  x  9  roosterpunten  voor 
i/8  kubus,  De  Boddingtonbenadering  levert  6dnzelfde  profiel  op  in  alle  richtingen, 
terwijl  dit  in  de  realiteit  en  ook  in  3D  ABAQUS  simulatie  niet  zo  is.  Van  de  ABAQUS 
simulatie  zijn  de  profielen  getekend  langs  de  as,  langs  de  01 1-diagonaal  en  langs  de 
1 1 1-diagonaal. 

De  Boddington  lijn  ligt  over  het  algemeen  een  paar  graden  boven  de  ABAQUS  lijnen. 
Bij  nadere  studie  blijkt  echter  dat  in  dit  geval  vergelijking  (A.  17)  voor  y„  geen  reele 
op  tossing  heeft  voor  n  =  2  en  oplossing  (A.  18)  is  in  dit  geval  dan  ook  niet  correct. 

Het  gevolg  is  dat  het  profiel  met  name  voor  p  — >  1  niet  klopt.  Hoewel  Boddington 
beweert  dat  de  2“  orde  benadering  al  goede  resultaten  geeft  blijkt  dat  voor  dit  specifieke 
geval  toch  niet  te  kloppen  en  is  een  hogere  orde  benadering  noodzakelijk. 


Tabel  A.  I  Temperatuurprofiel  kubus  voor  2'  nnle 


p 

0 

T(K) 

t  rej 

0,00 

1 ,8878 

367,03 

93,88 

0,04 

1.8814 

366,98 

93,83 

0,08 

1,8622 

366,82 

93,67 

0,12 

1,8308 

366,56 

93,41 

0,16 

1,7800 

366,20 

93,05 

0,20 

1t7348 

365,75 

92,60 

0,24 

1  f6723 

365,23 

92,08 

0,28 

1,6020 

364.64 

91.49 

0t32 

1,5251 

364.00 

90,85 

0,36 

1.4430 

363,31 

90,16 

0,40 

1,3571 

362,60 

89,45 

0.44 

1,2686 

361,86 

88,71 

!  0.48 

1,1786 

361,10 

87,95 

benadering  van  Boddington. 


P 

0 

T(K) 

T  (°C) 

0,52 

1,0882 

360,35 

87,20 

0,56 

0,9984 

359,60 

86,45 

0,60 

0,9100 

358,86 

85,71 

0,64 

0,8237 

358,14 

84,99 

0,68 

0,7403 

357,44 

84,29 

0.72 

0,6603 

356.77 

83,62 

0,76 

0,5844 

356,14 

82,99 

0,80 

0,5130 

355,54 

82.39 

0,84 

0,4466 

354,98 

81,83 

0,88 

0,3857 

354,47 

81,32 

0,92 

0,3307 

354,01 

80,86 

0,96 

0,2822 

353,61 

80,46 

1,00 

0.2407 

353,26 

80,11 
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n=2 


P 

Boddington  . Abaqus  as  (001) 

Abaqus  diagonaal  (01 1 ) - Abaqus  diagonaal  (111) 


Figutir  A.l  Temper  atuurprofielen  kuhus  vtmr  2e  t>rde  betiadering  van  Boddington  en  zoals  berckcnd  mol 
behulp  van  Abaqus- 


Bij  een  5U  orde  henadering  worden  de  kritische  waarden  gegeven  door: 

X(2n-l)6fl>;r  =  l  =>y«=  1,25625 

tt=  1 

=  2,0784 

n= 0 

0^  =  21n(xtJ  =  1.4632 
Ar  (^0  )  =  1  2  ~~Y~  =3-4898 

*  Ur 


TAtCr  =  350,82  K  =  77,67  °C 
Ttur  =  363,02  K  =  89,87  °C 

T^cjen  vooral  ToiCrzijn  nu  een  stuk  lager  dan  in  bet  geval  van  de  2C  orde  benadering. 
Het  temperatuurprofiel  voor  dit  geval  is  gegeven  in  label  A.l  en  figuur  A.l.  Hei  profiel 
ziet  er  nu  goed  uit,  houdt  ongeveer  het  midden  tussen  de  verschillende  door  ABAQUS 
berekende  proFielen  en  ligt  ongeveer  !  °C  onder  hel  door  ABAQUS  berekende  profiel 
langs  de  as. 
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Tabe  I  A .  2  T emperatu  u  ip  ro  fie  I  kubus  vtxi  r  5C  o  rde  be  naderi  ng  va  n  Bodd  ingn>n . 


p 

e 

T(K) 

TpC) 

0,52 

0*9065 

358.38 

85.23 

0*56 

0,8347 

357,78 

84*63 

0,60 

0,7615 

357,17 

84*02 

0,64 

0,6873 

356,55 

83,40 

0,68 

0,6124 

355,93 

82,78 

0,72 

0,5370 

355,30 

82.15 

0,76 

0,4614 

354*67 

81,52 

0,60 

0*3056 

354,04 

80,89 

0,84 

0*3096 

353,40 

80,25 

0*88 

0,2334 

352,77 

79,62 

0,92 

0,1566 

352,13 

78,98 

0.96 

0,0790 

351,48 

78.33 

1,00 

0,0000 

350,82 

77,67 

P 

0 

T(K) 

T  (°C) 

0*00 

1*4632 

363,02 

89,87 

0,04 

1,4592 

362,99 

89*84 

0,08 

1,4472 

362,89 

89*74 

0.12 

1,4275 

362,72 

89,57 

0,16 

1,4003 

362,50 

89,35 

0,20 

1,3660 

362,21 

89,06 

0,24 

1*3252 

361,87 

88,72 

0,28 

1,2784 

361,48 

88,33 

0*32 

1 ,2263 

361,05 

87*90 

0,36 

1,1694 

360,57 

87*42 

0,40 

1,1083 

360,06 

86*91 

0.44 

1,0438 

359,52 

86,37 

0,48 

0,9763 

358,96 

85,81 

n  =  5 


P 

Boddington  . Abaqus  as  (001) 

- Abaqus  diagonaal  (01 1 ) - Abaqus  diagonaal  (111) 


tfguur  A. 2  Temperatuurpmiieten  kubus  vnor  So  onio  bonadoring  van  Boddington  on  mala  horokond  mot 
hehuip  van  Abaqus* 
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